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ABSTRACT 

Watkins, William Allen. Ph . D . , Purdue University, 
December 1972. A Study of Lateral Diffusion in Thin Film 
Couples of Gold and Silver. Major Professor: 

Samuel J. Hruska. 

An experimental deposition technique was developed to 
model a thin film couple after the standard bulk diffusion 
sandwich couple. The initial region of mixing between the 
two sides of the thin film couple was less than 2 microns. 

An empirical conversion parameter, following the bulk 
procedures of Ziebold and Ogilvie, was determined for thin 

O 

films of gold and silver between 500 and 2,000 A. X-ray 
intensity data from the electron microprobe were converted, 
through use of the parameter, to units of mass concentration. 
Pure and alloy films of gold and silver exhibited linearity 
between the microprobe x-ray intensity and film thickness 

O 

up to 2,000 A. The maximum error in alloy film composition 
using the conversion parameter was shown to be 10 atomic % 
at 500 A. 

Thin film couples of pure silver and gold annealed above 
500°C, and at 416°C for times greater than 72 hours, exhibited 
dewetting on the pure gold side. The dewetting of the gold 
film was shown to be in accord with Mullin's theory of grain 
boundary grooving. For Ag-Au couples annealed below 500°C 
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;i zone of porosity occurred, increasing with time and tem- 
perature of anneal, leading to film separation near the pure 
silver side of the original region of mixing. The porosity 
formation was similar to that reported by investigators of 
bulk sandwich couples of pure Au and pure Ag . 

Incremental thin film couples with initial concentration 
differences of 20 atom% did not exhibit the zone of porosity 
which again was in agreement with investigators using bulk 
couples with small concentration differences. The diffused 
concentration profiles were compared with numerical solutions 
of the initial profile. The numerical technique was based 
on an approximate solution to Fick's second law for one 
dimensional diffusion. 

A series of apparent diffusion coefficients was deter- 
mined for temperatures between 200 and 416°C. A plot of 
In D versus 1/T was linear with a coefficient of correlation 
of 0.99. The data fit the following relationship 

D = 1.9 x 10" 7 exp (- 10 , 900/RT) cm 2 sec _1 

An apparent diffusion coefficient was also determined 
for a series of pure Ag and Au couples annealed at 416°C and 
compared favorably with the incremental couples' diffusion 
coefficient at the same temperature. 

Two models were considered in order to attribute the 
apparent diffusivity to a single transport mechanism. The 
surface diffusion model was refuted because of the low 
reported value of activation energy (10.9 kcal/mole), in 
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comparison with what was theoretically expected and experi- 
mentally reported by other investigators. A model attribut- 
ing the apparent diffusivity to transport along the film's 
grain boundaries was considered to be likely, if the grain 
boundaries could be assumed in a high state of disorder. 
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INTRODUCTION 

O 

Diffusion in thin (< 5,000 A) vacuum deposited metallic 
films is a phenomenon which has received significant atten- 
tion only during the last ten years. This is largely due 
to the increasing use of thin films in electronics either 
as complete thin film circuits, or more commonly as inter- 
connections and contacts in integrated circuits. 

Until the last two years, the majority of thin film 
diffusion studies had been done with multilayer films. 
Multilayer films are produced by the successive deposition 
of material stacked layer by layer, one on top of the other, 
as shown in Figure 1. In a layered film, deposited in high 
vacuum (< 10’ 5 torr), each interface region will be nearly 
perfect. For example, there will be atomic contact between 
any two superimposed deposits. Diffusion occurs in one 
direction only, and the limiting boundaries of the system 
are perpendicular to the direction of the highest concentra- 
tion gradient. Several investigators have devised experi- 
mental methods to detect the subsequent diffusion over the 
extremely short distances relative to that of bulk diffusion. 
Two comprehensive and fully referenced review articles have 
been written by C. Weaver [1, 2] and bring the reader up 
through 1970 in all work related to diffusion in thin films. 
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Figure 1 

Multilayer Film Geometry Using Alternating 
Films of Metal A and Metal B 
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Interconnections and contacts in integrated circuitry 
made of dissimilar materials have a possible mode of atomic 
transport which cannot be modeled using the multilayer film 

O 

technique. A thin film (< 2,000 A) couple made of two films 
side by side on a suitable substrate has an increase in the 
surface to volume ratio of approximately 10 5 relative to a 
standard size bulk diffusion couple. Atomic transport 
across the film surf ace - vapor interface could significantly 
contribute to the rate of homogenization, especially at 
temperatures below \ the melting point of the pure components 
of the couple. Film- substrate interactions and surface 
morphology changes for various annealing times and tempera- 
tures can also be studied. 

Lateral diffusion, diffusion along the length of a 
binary thin film couple parallel to the substrate, has not 
had the same degree of attention as that of diffusion normal 
to the substrate (multilayer films). Only four investigators 
[3, 4, 5 and 6] have developed experimental techniques to 
investigate lateral diffusion. Three authors [3, 4 and 5] 
used the overlapping film technique. This method consists 
of alternating film deposition onto a suitable substrate 
intervened by a mask shift parallel to the substrate. The 
resultant film-substrate geometry used by all four investi- 
gators is shown in Figure 2. 

After film deposition, each investigator gave the 
binary thin film couples an anneal in vacuo. The overlap 
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DIFFUSION DIRECTIONS 






A. 



B. 



C. 





REGION OF 
MIXING 



Figure 2 

Film Geometry for Previous Lateral Diffusion Studies; 

A. Monch [3] B. Pines et al . [4] and .Johnson and Brown [5] 

C. Smallen [6] 
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indicated by A. {j B. in Figure 2, homogenizes rapidly due to 
the short diffusion distance perpendicular to the substrate. 
Two of the four investigators [3, 4] overlapped films of 
Ag with Te and Sc, respectively. Another author [5] studied 
diffusion in both sets of binary systems. Both Ag-Te and 
Ag-Se have an intermetallic Ag 2 X 3 reported [7] in studies 
of their equilibrium phase system. All three authors [3, 4, 
5] reported rapid formation of the intermetallic compound in 
the overlap region. Further, each of these investigators 
studied the subsequent lateral growth of the intermetallic 
after homogenization of the overlap region. High rates of 
diffusion with low activation energies, relative to bulk 
diffusion values, were reported. [3, 4] attributed these 
results to transport along grain boundaries, and [5] to 
diffusion along the film surface- vapor interface. In all 
these three independent studies, the experimental results 
were strongly dependent on the subsequent film geometry. 
Variations in either of both of the film's thicknesses 
caused variations of two orders of magnitude in the reported 
diffusion rates. In [5], the Ag had to be deposited first 
on the substrate in order for lateral diffusion to occur at 
all . 

Smallen's work [6] in the gold-copper system provided 
the major impetus for the present study. Using the method 
of simultaneous deposition of the gold and copper onto a 
suitable substrate, he reported rates of diffusion both 
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parallel and perpendicular to the substrate. Smallen 
reported rates of diffusion in the lateral direction almost 
seven orders of magnitude faster than in the perpendicular 
direction. An associated activation energy was reported and 
was almost | that of the previously reported bulk value. 

A study of the Bi-Te system, which also has a reported 
[8] intermetallic compound in its equilibrium phase system, 
was undertaken by this author. In order to more closely 
approximate the sandwich couple, so as to compare more 
favorably with the initial conditions of the bulk, it was 
decided to try to improve the experimental methods of [3, 4, 
5 and 6] . 

An experimental method was devised and reported [9] in 
order to reduce the overlap region of [3, 4, 5 and 6] from 
the order of 10' 2 meters, to 10' 5 meters, and minimize the 
thickness variation in the lateral direction. The new 
geometry for a thin film diffusion couple used in the 
author's work is shown in Figure 3. Note that the overlap 
region of Figure 2 is now called a region of mixing, which 
more closely represents the physics of the new experimental 
deposition process [9] . 

The formation of an intermetallic compound in the 
region of mixing was observed in couples made of pure 
bismuth and tellurium. However, after formation of the 
intermetallic compound, no lateral diffusion was observed. 

A summary of the work accomplished on this system, with 
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REGION OF MIXING 
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DIFFUSION 

DIRECTIONS 



Figure 3 

Geometry for Thin Film Couples Used in Present Study 



the experimental results, and a discussion of the same, is 
found in Appendix A. 
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A second study was undertaken in the Ag-Au system, 
which is entirely miscible throughout the entire composition 
range. In this study, the region of mixing was further 
reduced to the order of 10" 6 meters as the experimental 
procedure was further improved. (See Experimental Procedure 
section . ) 

The initial phase of the second study involved couples 
of pure Au and Ag . During the anneals surface morphology 
changes occurred in the couple. The changes were similar 
to those previously reported in studies of bulk Ag-Au 
couples, however they occurred to a greater extent and 
interrupted the continuity of the couple. Because of the 
film separation, quantitative rates of diffusion could not 
be determined as anticipated. However, similarities between 
the concentration profiles for various annealing times led 
to an alternate technique; this allowed values of diffusivity 
to be determined. 

Johnson [38] , using bulk Ag-Au couples with small 
initial concentration differences, did not observe these 
morphology changes. Accordingly, an experimental procedure 
was devised to make thin couples with an initial concentra- 
tion difference of 20 atomic %. As with Johnson's results, 
no morphology changes were observed in the present study. 
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A numerical technique was devised to treat experi- 
mentally determined initial concentration profiles. The 
numerical solutions were compared with the experimental 
concentration profiles after various times and temperatures 
of anneal. Values of chemical diffusivity were determined 
as a function of annealing temperature. From the comparison 
of the numerical solutions with the diffused concentration 
profiles, an equation for the apparent diffusivity for the 
couple was determined. Possible atomic transport mechanisms 
were investigated over the working range of the present 
s tudy . 

Lastly, a calibration study, involving the fabrication 
and use of Ag-Au thin film standards was undertaken. 

Lmpirical working curves over the entire composition 
range were determined so that electron micro-analysis could 
be directly applied to deduce mass concentration variations 
in the diffusion studies. 
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BACKGROUND AND LITERATURE SURVEY 

Two separate investigations into the literature had 
to be undertaken in order to carry out the experimental 
work of this author. 

The first investigation focused on a method to analyze 
the annealed thin film couple and from this, deduce mass 
concentration variations as a function of lateral distance, 
time, and temperature of anneal. 

The second involved an investigation into the previous 
diffusion studies in the Ag-Au system in order to compare 
with the experimental results of the present thin film study 

Empirical Me thod for Thin Film Analysis 

The experimental method of [3, 4, 5] involved the 
optical observation of the growth of the diffusion zone as 
a function of time and temperature. In all three cases, the 
diffusion zone was an intermetall ic compound with a signifi- 
cant reflectivity difference relative to the pure components 
Smallen [6], however, worked in the gold-copper system which 
is completely miscible throughout the composition range. 

In order to determine atomic concentration profiles in his 
couples he used the electron microprobe. By scanning his 
samples, under the electron beam, he generated x-ray 
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intensity data as a function of lateral distance. In order 
to convert the x-ray spcctromc tr i c data to atomic concen- 
tration, he made the rather crude assumption that the 
relative x-ray intensity, corrected for background, was 
proportional to atomic concentration. 

Due to the similar work in miscible systems by Smallen 
[6], and to the present work of this author, the electron 
microprobe was used to determine mass concentration changes 
as a function of lateral distance, time, and temperature of 
anneal. However, an investigation was necessary to deter- 
mine the relationships between the x-ray spectrometric data 
generated from a thin film sample, under the electron beam, 
and the specimen's mass concentration. 

The important variables relating the emitting x-ray 
intensity to the chemical composition of the beam target, 
were first reported by Castaing [17]. The full correction 
for the measured relative intensity in the bulk is 



K 



A 



^alloy 

I A 

pure A 



f(x) 

Lf°(x)J 



1 



+ 




"RE_J 

_R°p°J 



J A 



CD 



where [f(x)/f°(x)] 
[1 + K £ ] 



is the absorption correction factor 
is the correction for secondary 
fluorescence 



R P 

lrVj 



is the backscatter- ionization correction 



C 



A 



is the weight fraction A in A-B alloy 



12 



lalloy * s t * ie measure d x-ray intensity of the 
characteristic line for element A 
(corrected for background) of the A-B 
alloy 

Ipure a t ^ ie measure d x-ray intensity of the 

characteristic line for element A 
(corrected for background) of a sample 
of pure A. 

Ziebold and Ogilvie [18] made a series of experimental 
observations between and in 17 different binary 
systems, including the Ag-Au system. They found that their 
experimental data in all 17 systems displayed linearity 
when plotted in the form of C^/K^, versus C^. The analyti- 
cal straight line function is written 



c a /k a 



= a 



AB 



(1 - 



Vb^a* 



( 2 ) 



The single parameter needed to describe the empirical work- 
ing curves of [18] is the limiting value of C^/K^, when 
goes to zero. The subscript, AB , indicates that the con- 
version parameter applies to the measurement of characteris- 
tic x-rays of element A in binary A-B alloys. 

Rearrangement of equation (2) gives a more symmetrical 

form: 



1-K 



A 



1-C 



K 



= a 



A 



A 



AB 



(3) 
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Figures 4 and 5 show the experimentally determined bulk 
relationships of [18] between K and C in the silver-gold 
system. A hand calculation of the values of the three 
correction parameters of equation (1) , as a function of 
composition, was made using semi- theoretical formulae, 
based on the three correction factors of equation (1) . 

An operating voltage of 25 kv of the electron microprobe, 
with a 52.5° spectrometer takeoff angle, (parameters used 
in the present experimental work) was assumed in the calcu- 
lation. In Figure 6, the results of the hand calculation 
are plotted in terms of the variables of equation (3) . 

In the case of thin film analysis, if the linear 
analytical function of equation (2) is assumed to hold 
for thin films in the Ag-Au system, then the empirical con- 
version parameters (a^gAu’ a AuAg^ must be determined 
separately from the previously reported bulk values. 

Hutchins [19] showed that the thin film x-ray intensity, 
normalized with respect to a bulk sample of the film material, 
can be expressed as three factors: 



= (S) (i F ) 

DUlk 



1 



I 



bulk exc. 






(4) 



where Ipg = the measured intensity for the film on a 
particular substrate 

*bulk = the measured intensity for an infinitely thick 
film (with respect to the electron beam) 

S = substrate correction 



MEASURED RELATIVE INTENSITY 
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Figure 4 

Experimental Microprobe Calibration Data for the Bulk 
Silver-Gold System at 30 kv and 15.5° Spectrometer 

Takeoff Angle. 

Zeibold and Ogilvie [18] 
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(Ca) 

Figure 5 

Functions C^/K^ for Bulk Silver-Gold System. 
Microprobe at 30 kv and 15.5° Spectrometer 
Takeoff Angle. 

Zeibold and Ogilvie [18] 
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Figure 6 

Results of Hand Calculation Using Semi-Theoretical Formulae, 
Based on Equation (1) Plotted in Form of Equation (3) . 
Microprobe at 25 kv with 52.5° Spectrometer Takeoff Angle. 
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Ip = the measured intensity for the film without a 
substrate 

I bu 1 k exc = the excited intensity for an infinitely 
thick film 

f(x) = the x-ray absorption correction for film 
thickness, t. 

For a particular set of operating conditions of the 
electron probe and film thickness, t, the measured relative 
intensity can be written as 



1 alloy(FS) 
Ipure A(FS) 



(5) 



where I a ]_]_ 0 y(pg) = the measured intensity for an A-B alloy 

film on a particular substrate at a film 
thickness t 

Ipure A(FS) = t * ie measure d intensity for a film of 
pure A on a particular substrate at a 
film thickness t. 

Solving equation (4) for Ipg, and substituting into the 
numerator and denominator of equation (5) with appropriate 
subscripts 



K^(film thickness 

Ipure A bulk exc. 
J^alloy bulk exc. 



t) 



I alloy (F) . 
Ipure A(F) 



f(x) p 

f (x) A 



*bulk alloy 
^bulk pure A 



( 6 ) 
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Ziebold and Ogilvie [18] were unable to justify their 
experimental observations leading to equation (2) on the 
theoretical grounds of the bulk relationship of equation (1) . 
This author will likewise make no claim to the prediction of 
a linear relationship in the present thin film data plotted 
in the form of equation (2) , based on equation (6) . A 
straight line relationship, within the limits of the repro- 
ducibility of the experimental data, between the variables 
of equations (2) or (3) will justify the use of the empirical 
conversion parameter to deduce weight fraction composition 
from the measured relative intensity of the electron micro- 
probe . 

To convert from weight fraction composition (C^) to 
mass concentration involves knowledge of the density varia- 
tion of gold-silver alloys across the entire composition 
range . 

In Figure 7, the Vegard's law deviation is shown for 
the silver-gold system. Using the lattice parameter data of 
Sachs and Werts [20] , and the relationship derived in 
Barrett and Massalski [22] for the density in a substitutional 
alloy system, 

p(AgAu) = nA/VN (7) 

where N = Avagadro's number 

V = volume of the unit cell 
N = number of atoms/unit cell 
A = mean atomic weight, 



LATTICE PARAMETER (A) 
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Figure 7 

Variation in Lattice Parameter of Bulk Alloys in the 
Silver-Gold System Using Experimental Data of 
Hund and Trager [21] and Sachs and Weerts [20], 



the density of gold-silver alloys as a function of atomic 
concentration was determined and is shown in Figure 8. 
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Bulk Diffusion in the Gold - Silver System 
At the present time there has been no reported study of 
thin film diffusion in the gold-silver system. 

The initial studies of chemical bulk diffusion in the 
silver-gold system was undertaken by [23, 24 and 25]. In 
all three works, the method of chemical etch limit was used 
to determine average diffusion coefficients over the entire 
composition range. Johnson [38], who doubted the validity 
of the chemical etch limit method, studied chemical diffusion 
in the silver-gold system using sandwich couples with initial 
concentration differences of less than 20 atomic % . Concen- 
tration-penetration profiles were determined by the now 
accepted method of sectioning. Using Fick's Second Law, 



3C = 3_ f n ^C' 

9t 9x 9x - 
\ / 



( 8 ) 



Johnson assumed that, due to the use of couples with small 
initial-concentration gradients, D was not a function of 
composition over the working range of his couples . Therefore 
he wrote 



9C _ n 9 2 C 
at " " 9x 2 



( 9 ) 



and solved equation (8) using the standard initial and 
boundary conditions. 
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Figure 8 

Variation of Alloy Density for Bulk Alloys in Gold-Silver 
System Using Relationship Derived in Barret and Massalski [22] 
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C -C 0 = 1 
^ 1 " ^ o 2 



1 - erf 



2 ^ 



( 10 ) 



where C = concentration, after time t, at a distance x 
from the interface 

C x = concentration at t = 0 , for x < 0 
C 0 = concentration at t = 0, for x > 0 



Johnson compared his experimental profiles with those of 
equation (10) to determine D, and also plotted his results 
on probability paper to justify his assumption of D f f(C). 
Figure 9 shows Johnson's results for an alloy containing 
49.2 atomic I gold, along with the results of [23, 24 and 
25]. Jost's [25] data, taken at lower temperatures, were 
extrapolated by Johnson [38] up to his experimental tempera- 
ture range. The data of an investigator named Braune are 
reported in Johnson's work and shown in Figure 9, however, 
he does not reference him in his paper. The straight line 
through Johnson's data can be represented by the following 
equation , 



D 

c 



0.14 



exp 



[-41,700] 

RT 

V / 



cm 2 sec' 1 . 



(ID 



Included in his work were determinations of the gold 
and silver self diffusion coefficients using radioactive 
tracers. Over the same temperature range they are repre- 
sented by the following equations. 



IN CM /SEC 



°C 




1000 

~ T“ 



Figure 9 

Chemical Diffusivity as a Function of Temperature 
for Bulk Sandwich Couples in the Silver-Gold System. 

Johnson [38] 
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D Au = 
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exp 


[-44,100] 
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(13) 



A very extensive study of self diffusion in silver-gold 
solid solutions was done by Mallard et al . [26], also using 

radioactive tracer sectioning techniques . They reported 
that the silver tracer moved between 2 and 3 times faster 
than the gold tracer, across the entire composition range. 
They also report a variation of less than \ of an order of 
magnitude in the dif fusivities of both the gold and silver 
tracers across the entire composition range. 

Baluffi and Seigle [27] studied diffusion in bi-metal 
vapor-solid couples of zinc-alpha brass, copper-nickel and 
gold-silver. In their investigations of vapor-solid couples 
of gold-silver they observed the formation of 13% volume 
porosity, by lineal analysis, in the diffusion zone after 
93 hours at 940°C. In order to elucidate the cause of the 
porosity formation, they made an additional study of pure 
gold-pure silver sandwich couples. The porosity again 
appeared in the sandwich couple and its formation was 
evident at times as early as 1 hour, increasing with time 
up to 50 volume %. Similar observations on gold-silver 
sandwich couples were made by Seitz and Kotterman [28]. 

Both investigators [27, 28] observed the porosity formation 
on the pure silver side of the couple. Baluffi and Seigle 
[27] made corrections to account for the observed volume 
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porosity to determine dif fus ivities by the Matano method. 
Seitz and Kotterman did not make corrections and therefore, 
their data are suspect. Baluffi and Seigle took Laue back 
reflection x-ray photograms of the diffusion zone. They 
observed multiple spots and attributed this to polygoniza- 
tion of the diffusion zone. 

In his earlier work, Baluffi [29] determined that 
silver atoms diffuse into gold more rapidly than gold 
diffuses out, causing a volume increase, with expansion in 
all three directions. Using the theoretical calculations 
of Huntington and Seitz [31] , he speculated that diffusion 
by vacancy mechanism coupled with the unequal diffusion 
rates, could account for the observed volume increase. 

Barnes [30] used multiple sandwich couples of thin 
interleaved sheets of copper and nickel to accentuate any 
observed volume increase due to diffusion. He reported that 
there was an increase of volume in the diffusion zone as a 
result of atomic diffusion in the copper-nickel and copper- 
zinc systems that he studied, each system having unequal 
rates of the two diffusing species. He also showed that 
this volume increase is associated with the formation in 
the diffusion zone of voids with faceting. These voids 
appear near the original interface and on that side from 
which there is a net loss of atoms (net gain of vacancies) . 
He also conducted an experiment which showed the presence of 
strain in the diffusion zone and another experiment to 
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attribute this strain to diffusion. He conducted an x-ray 
investigation like that of Baluffi and Seigle [27], and 
again concluded the presence of polygoni zation in the diffu- 
sion zone. Barnes [30] speculated that the difference in 
atomic fluxes of the two diffusing species was compensated 
by a vacancy flux in the direction of the smaller flux. 

The generation and annihilation of vacancies on opposite 
sides of the couple cause strain in the couple. He sug- 
gested that the voids forming as a result of diffusion, 
imply that some of the large number of vacancies involved 
in the diffusion zone condense out of the lattice. He 
postulated generation and annihilation sites such as grain 
boundaries, polygon boundaries and on the edge component of 
dislocations . 

He, however, did not discuss the magnitude of the devia- 
tion of the equilibrium number of vacancies present in the 
diffusion zone during the anneal. He explained that the 
growth of the void could be explained by diffusion and 
absorption of some of the excess numbers of vacancies on 
void surfaces . 

Seitz [34] postulated that the process of porosity 
formation is precipitation of supersaturated vacancies 
caused by the unequal diffusion current and occurs by 
heterogeneous nucleation. He calculated a value of 10 11 for 
the lifetime number of vacancy jumps in the diffusion zone 
by using the porosity observations of Barnes and an arbitrary 
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value of 1 for the critical relative excess vacancy concen- 
tration to form holes. 

Baluffi [32] supported Seitz's postulate of heterogene- 
ous nucleation to form voids from a supersaturated solution. 
Baluffi, however, pointed out that the high value of super- 
saturation expected by Seitz could not be compatible with 
the experimental results of parabolic Kirkendall marker 
shifts with time. He considered the problem with an 
improved mathematical treatment, and his results indicated 
that the relative excess vacancy concentration for void 
formation to be < 0.01 (II). 

Barnes and Mazey [33] using sandwich couples of Cu and 
Ni under a hydrostatic pressure reaffirmed Baluffi 's [32] 
conclusion of II relative excess vacancy concentration 
necessary to form holes, and also showed experimental evi- 
dence that void formation was by heterogeneous nucleation, 
probably foreign particles or inclusions at the weld inter- 
face. A more recent and complete mathematical treatment of 
the vacancy flux contribution during diffusion and Kirkendall 
shift in binary alloys, is given by Manning [37]. 

During an experimental investigation of the Kirkendall 
shift in the gold-silver system, Ruth [35] observed the 
formation of a zone of small protrusions of the pure gold 
side of the couple. He reported that the growth law of 
the width of this zone is a secondary effect of diffusion. 

The protrusions were observed to form in the vicinity of 
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intersections with crystal surfaces of <112> lines in {110} 
planes. He explained the observations by proposing that a 



acting as a diffusion pipe. The excess silver atoms, 
caused by the unequal diffusion rates, should migrate along 
the dislocation line to the surface. 

There are no reported studies of grain boundary diffu- 
sion in bulk silver-gold alloys. Hoffman and Turnbull [36], 
however, studied lattice and grain boundary self diffusion 
in silver. Over the temperature range of their study (250- 
500°C) they reported the grain boundary diffusion coeffici- 
ent times the grain boundary thickness fit the equation, 



^ [110] edge dislocation, which has a [112] direction, 




where y = 5 * 10 ‘ 8 cm 



Grain sizes for their specimens varied from 21-60 microns. 
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DESIGN AND CONSTRUCTION OF EQUIPMENT 

The successful deposition of thin films in vacuum 
requires a working understanding of vacuum technology. The 
experimental apparatus must be maintained in an environment 
suitable for the selective transport of mass over a distance, 
using materials which maintain this environment during all 
the operating conditions of the experiment. In order to 
give the reader insight into the many criteria which had to 
be satisfied in order to successfully create the three thin 
film geometries, and carry out the required experimental 
conditions for the study of atomic diffusion, detailed 
explanation is given in the following sections. 

The discussion of this author's experimental equipment 
is broken into three sections: the vacuum system, the 

experimental apparatus in vacuo, and the associated external 
equipment . 



Vacuum System 

The environment for the generation of both binary thin 
film couples and thin film standards was a vacuum chamber 
capable of maintaining residual gas pressures of 10‘ 7 torr. 
Primary pumping was obtained using a 4 inch oil diffusion 
pump in conjunction with a mechanical roughing pump. 
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Auxiliary pumping, located in a stainless steel base 
separated from the diffusion pump by a high vacuum valve, 
included a Varian tri-filament titanium sublimation pump 
and a liquid nitrogen cooled cyrogenic baffle. Residual 
gas pressure was measured using a Granville-Phillips dual 
tungsten filament ionization gauge tube with an associated 
Series 260 controller. 

A 12" Varian, all stainless steel, spool piece was 
placed on top of the base, using viton gaskets to maintain 
vacuum integrity. Integral to the spool piece are ten 2jj" 
Conflat ports for insertion of necessary feedthroughs to 
maintain electrical and mechanical connectivity between the 
evacuated apparatus and the external test and monitoring 
equipment. An 18" high pyrex glass bell jar, also sealed 
with a viton gasket, was placed on top of the spool piece 
completing the vacuum system described above. A sketch of 
the vacuum system with its important components labeled is 
shown in Figure 10. Photographs of the exterior and 
interior of the vacuum system are displayed in Figures 11 
and 12 respectively. 

Experimental Apparatus in vacuo 

The substrate holder, two sources, and associated monitor- 
ing and shielding mechanisms were constructed entirely within 
the 12" diameter spool piece. OFHC copper, type 303 and 304 
stainless steel, aluminum, and 96% polycrystalline dense 
alumina, were used exclusively in the construction of the 
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Figure 10 

Schematic of the Vacuum System Used in the Present Study 
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Figure 11 

Photograph of the Vacuum System and Associated 
External Equipment Used in the Present Study 
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Figure 12 

Photograph of the Experimental Apparatus in vacuo 
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experimental apparatus. These materials have reported [10] 
low (< 10' 9 torr) vapor pressures at the operating tempera- 
tures of the vacuum system (-192°C to 600°C) . In addition, 
wherever possible, leak paths were provided for all fastened 
assemblies to minimize the outgassing of virtual leaks. 

The substrate holder, a sketch of which is shown in 
Figure 13, was made from a 3" x 2" x I" block of OFHC copper. 
Both large faces of the copper block were milled to within 
.001". Two similarly milled copper restraining plates were 
used to hold the substrate uniformly against the copper 
block. All materials, except the heater and insulating tubes, 
which will be described later, were constructed out of OFHC 
copper. The use of copper exclusively on the substrate 
holder eliminated differential stress buildup during heating 
of the block, maintaining the critical alignment between the 
various support mechanisms both on and off the substrate 
holder . 

A copper quench tube, making two passes through the 
substrate holder, was connected to two stainless steel Swage- 
lock fittings. The fittings completed the quench line, 
through attachment to two stainless steel quench gas delivery 
tubes connected externally to the gas source by a stainless 
steel dual tube feedthrough. In addition, the fittings 
allowed removal of the substrate holder, by sealing together 
through the use of replaceable copper gaskets. The razor 
blade support mechanism was made of two sets of two small 
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Figure 13 

Schematic of Bottom View of the Substrate Holder 
with Substrate and Razor Blade Mounted in Place 




37 



copper blocks, squared true to the substrate holder's milled 
face. In each set where the faces held the razor blade, a 
.002" reset was milled on each face to accommodate the razor 
blade thickness. Copper screws held together the two sets 
of blocks in each set, and as a unit to the substrate holder. 

Chromel - alumel thermocouples used to record the sub- 
strate temperature were press fit by screw heads onto the 
surface of the substrate away from the line of sight of the 
two sources. 

Four (4) tubes, symmetrically placed above and through 
the substrate and copper block respectively, were filled 
with sleeves of 99.7 wt. % alumina of 96% polycrystalline 
density. The sleeves provided electrical insulation between 
the copper block and the .010" coiled tungsten wire, which 
was interwoven between the four sleeves, providing a heat 
source to the substrate. Alumina of lesser 
density was found to absorb water vapor readily at atmos- 
pheric pressure and subsequently outgas heavily during 
heatup in vacuum. 

The deposition of gold and silver was carried out in 
two independently controlled resistance heated boats 
purchased from the R. D. Mathis Company. A .010" molybdenum 
boat with a .010" coating of alumina was used for the gold 
deposition, and a .010" tantalum boat was used for the silver. 
The tungsten boat more commonly used for gold deposition in 
flash evaporation systems was found to be partially wet by 
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the gold, giving erratic rates of deposition, and possible 
stratification of the deposited film. 

Figure 14 is a sketch of the approximate orientation 
relationships between the sources, substrate holder, shutter 
mechanism, and crystal mass deposition monitors. 

The shutter mechanism moves, according to the plane of 
the sketch in and out of the paper, on two aluminum tracks 
fastened to the spool piece. The shutter was designed with 
the winged contour to provide, simultaneously, controlled 
interruption of the gold and silver vapor to the substrate, 
continuous monitoring of each vapor flux by the respective 
crystal monitors, and vibrationless transit with respect to 
the substrate-razor blade-flux barrier alignment. 

As indicated in Figure 14, the flux barrier between the 
two sources is completed by a stainless steel sheet. In 
order to compensate for expansion of the substrate holder 
during heating, and adjustment of the edge of the sheet true 
to the edge of the razor blade, the flux barrier is adjust- 
able vertically and horizontally, respectively. 

The radiant transfer of heat to the vacuum chamber 
walls and surrounding apparatus during heating of the sub- 
strate holder assembly was minimized by the construction 
of a two piece (for disassembly) copper box surrounding the 
substrate holder on three sides (open side facing sources) . 
The addition of the copper shroud served a dual purpose. 

It maintained the internal surfaces of the vacuum surfaces 
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Sketch of the Approximate Orientation Relationships 
between the Sources, Substrate Holder, Shutter Mechanism 



and Crystal Mass Deposition Monitors 
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below 100°C, preventing outgassing of the same reducing the 
low pressure capability of the system. Secondly, it assisted 
in maintaining an isothermal condition for the substrate. 

Two Sloan water cooled ultra-high vacuum remote crystal 
holders were rigidly mounted in line of sight, one to each 
source. Five (5) MC gold and silver coated quartz, AT cut 
round crystals were used in the holder for monitoring the 
gold and silver source respectively. Each crystal gives 

O 

thickness monitoring up to 60,000 A of aluminum or mass 
equivalent. The round crystals used were found to be 
superior to the square crystals used in early experimenta- 
tion. The uniform stress pattern generated in the round 
crystal during mass deposition eliminated false (unrelated 
to mass deposited) frequency shifts common to the square 
type. A radiant heat shield placed around the crystal 
minimized false frequency shifts of the crystals due to 
induced thermal stress. 

External Equipment 

Automatic deposition control was achieved using two 
Sloan Omni II Deposition Control Monitors. While each 
source was controlled to give the preset deposition rate, 
the rate and thickness output was continuously monitored 
by 4 Rustrak chart recorders. 

Substrate temperature was monitored by a Leeds and 
Northrup millivolt potentiometer, Sargent Model SRG chart 
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recorder and a Leeds and Northrup M Line Model C-l, C.A.T. 
controller. 

In conjunction with the C.A.T. controller, a Leeds and 
Northrup M Line Model C-2 low level set point unit, silicon 
controlled rectifier, and tungsten heating element in the 
substrate holder assembly was used to control the substrate 
temperature during the anneals to within + 1°C. 

A variac auto transformer fed by 115 VAC from a con- 
stant voltage transformer was used for initial heatup and 
outgassing of the substrate and was switched in or out of 
the heating circuit with a DPDT high current switch. 

The bottled helium quench gas was of commercial purity, 
regulated, and provided to the appropriate feedthrough, 
maintaining the substrate holder assembly, by tygon tubing. 
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SUBSTRATE 



When studying the physical phenomena in thin films, 
the substrate or restraining medium on which the films are 
created is of signal importance. The substrate in many 
cases may determine whether or not a particular physical 
phenomena in the film or films occurs. It was with this 
in mind that this author has written this section, and 
divided it into two parts, the selection of the substrate, 
and the subsequent substrate preparation prior to actual 
film deposition. 



Substrate Selection 

Single crystal silicon wafers grown in the [111] from 
the melt were purchased from the General Diode Corporation. 
The wafers were 1.5" in diameter and 10-15 mills in thick- 
ness. As received, one side was mechanically polished, the 
other etched. 

The selection of the silicon wafer as a substrate for 
all experimental work in the silver-gold system was based 
on the following criteria: Firstly, the substrate must 

have a surface flatness and roughness so to consistently 
reproduce the desired film geometry with respect to the 
region of mixing and uniform thickness. Deviations from 



43 



flatness in the substrate leads to variations in the distance 
between the razor blade edge and substrate. This coupled 
with the intrinsic variations in the razor blade edge height 
can lead to umbra effects resulting in formation of a gap 
in lieu of a region of mixing between the two films. The 
surface roughness of a substrate leads to surface protrusions 
of high frequency. The subsequent decay by surface self 
diffusion is considered an undesirable variable for the 
present experimental method. Doctors Anderson and Neudeck 
[11] compared nine common substrate materials and found the 
mechanically polished silicon wafer to exhibit a flatness 
on the order of 0.1 part/ thousand , and a surface roughness 
on the order of, or below, the resolution limit (* 50 X) of 
the measuring instrument [12] . 

Secondly, the thickness of the substrate was to be 
minimized as much as possible to reduce the temperature 
gradient maintained between the substrate holder where the 
temperature was measured and the surface of the substrate 
supporting the films. A 1.6 mm thick soft glass substrate 
clamped to the smooth face of a copper block held at 90°C, 
exhibited a 20°C temperature difference across the substrate 
thickness [13]. The silicon wafer, with 1,500 X of SiC >2 
grown on the surface, is almost an order of magnitude 
thinner than soft microscope glass. A temperature difference 
of 1°C, or less, across the wafer's thickness was experi- 
mentally observed by this author when the substrate holder 
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was held at 416°C, and uniform clamping pressure by the Cu 
restraining plates was applied. 

Thirdly, a substrate with a low average atomic number 
was desirous for quantitative thin film measurement by 
electron probe microanalysis. High energy electrons back- 
scattered from the substrate and the continuous spectrum of 
the substrate, will enhance the characteristic (Au and Ag) , 
x-ray intensity produced by the primary electron beam. Both 
the backscattered electron fraction and continuous spectrum 
intensity increase with atomic number [15] . The enhancement 
of the background intensity limits the theoretical limit of 
detection of trace elements. Campbell [14], using fluores- 
cent x-ray spectrographic analysis defined the minimum 
detectable limit as "that concentration or amount required 
to give a line intensity above background for counting times 
of 10 minutes." For smaller counting times, since the 
standard deviation is a square root function, the loss in 
sensitivity goes as a square of the factor decrease in count- 
ing time. Birks [16], comparing the methods of Campbell 
with the electron probe, showed that due to the smaller 
counting rates of the probe the limits of detectability are 
2 orders of magnitude higher (1 ppm vs. 100 ppm). 

An additional consideration for electron microanalysis 
of thin films is the presence of materials in the substrate 
which cause fluorescence excitation by line radiation. Tin 
and thorium generate a high fluorescent yield for Ag and Au 
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respectively, however these materials are not found as major 
constituents in most common substrate materials. 

fourthly, the linear coefficient of expansion of the 
substrate must be as close as possible to those of Ag and 
Au to prevent shivering or crazing to develop in the film, 
during the quench to room temperature at the conclusion of 
the anneal. Silicon, while not an ideal match, with a 
linear coefficient of expansion of 7.3 microinch/inch°C , 
vice 14.2 and 19.7 for the gold and silver respectively, 
gave no experimental observations of shivering or crazing 
for quenches up to 500°C. 

The above criteria are all satisfied by the selection 
of a silicon crystal wafer as a substrate. Individual 
criteria may be met to a greater degree by a different type 
of substrate, however none was found that satisfied all four 
criteria to the degree as that of the silicon wafer. 

Two criteria were not met by the silicon wafer; however, 
they were satisfied by the growth of a thin film of SiO^* 
by means of oxidation of the wafer. 

The growth of the amorphous layer of SiO^ allowed the 
deposition of polycrystalline thin films. Secondly, it acts 
as a diffusion barrier, allowing negligible diffusion perpen- 
dicular to the substrate for the times and temperatures of 
this work. 

* The oxidized layer contains molecules of SiO, SiC> 2 , and 
loosely bonded oxygen radicals. 
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Substrate Preparation 

All substrates were cleaned by immersion in hot 
sulfuric and nitric acid for 15 and 10 minutes respectively. 

A rinse in deionized water was followed by several rinses 
in reagent grade acetone and trichloroethylene. The sub- 
strates were then placed into a constant temperature zone 
of a tubular furnace fed with dry oxygen at 1,000°C, where 

O 

1,500 A of SiO were grown onto the silicon wafer. From the 
sensitive color change due to selective light absorption of 
the oxidized layer, its thickness was visually determined to 

O 

be uniform over the wafer to well within 100 A. 

The wafers were then masked and placed into a MRC 

O 

sputtering chamber, where a 5,000 A thick annulus of tantalum 
with an inner and outer diameter of 1" and 1.5" respectively, 
was sputtered onto the oxidized silicon wafer. The sputter- 
ing parameters were as follows: 



System Base Pressure 
Sputtering Gas 
Sputtering Gas Pressure 

Cathode Voltage 
Current 

Deposition Rate 



10" 5 torr 
Argon 

8 x 10‘ 2 torr 
(pressure controlled) 

4 kv D.C. (regulated) 

90 milliamps 

250 A/minute 
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The 5,000 A thick annulus maintained a fixed distance 
between the razor blade edge resting on the tantalum surface 
and the SiO surface. This fixed distance for each experi- 
mental run gave a reproducible region of mixing in the 
resultant film geometry. Figure 15 shows the top and edge 
view of the substrate after the preparation process 
described above. 



48 



A. 



B. 






BEVELED 
EDGE “ 

Si O 2 



RAZOR BLADE 



YZ TZZ ZZ ZZ Z ZZ TZZ ZZZZ ZZ Z Z 



To 



I THIN FILM j 



To 



SILICON WAFER 

[111J 



1.5 INCH 



5,000 & 
1,500 & 



T 

10-15 MILLS 

_i 



Figure 15 

Orientation Relationships between Silicon Wafer and 
Associated Films. A. Top View of Substrate 
B. Cross Sectional View of Substrate with 



Razor Blade in Place 
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EXPERIMENTAL PROCEDURE 

Three different film deposition methods were used during 
the course of the experimental work of this author. The pro- 
cedures used for the vacuum deposition of thin film standards, 
thin film couples, of pure gold and silver, and an incremental 
couple with an initial concentration difference of 20 atomic 
%, will be described in three subsections. 

Thin Film Standards 

Oxidized Silicon wafers, with no Ta annulus sputtered 
on the surface, were used as substrates. The 1.5" wafers 
were cut into small sections using a diamond scribe, and 
mounted onto the substrate holder. The available area of 
deposition was further reduced by masking with a thin .010" 

Ta sheet. This provided a sharp edge for thickness measure- 
ment, and a negligible variation in thickness over the film 
area. The razor blade was not used and the stainless steel 
barrier was lowered away from the substrate holder, its only 
use was to provide a line of sight barrier between the two 
sources to prevent cross contamination. The system was then 
evacuated to 5 x 10 -7 torr while the substrate was outgassed 
by heating to 150°C, driving off adsorbed water and other 
volatile impurities on the surface. 
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The sources were heated slowly over a period of 5 
minutes, followed by a 2 minute soak cycle just below 
deposition temperature to bring each source to a uniform 
temperature prior to deposition. 

Deposition was carried out simultaneously, the rates 

O 



(A/sec) of each source adjusted, prior to opening the 
shutter, by the control monitors, so that the time of 
deposition for both sources was equal, and yet still main- 
taining the necessary mass ratio to give the film concentra- 

O 

tion desired. Films of 500, 1,000, 1,500 and 2,000 A in 
thickness, each with compositions of 0, 20, 40, 60 and 100 
atomic %, were deposited. A detailed explanation of the 
method of achieving film composition is given in Appendix B. 

Gold and silver of 99.99% purity were used as source 
material for deposition of thin film standards. 

Film thickness was measured by a Dektak. Four (4) 
measurements of the edge height were taken on each film, 
with the average thickness reported. 

X-ray emission data were taken with the ARL-AMX Electron 
Microprobe. Intensity measurements were obtained with the 
electron microprobe, operated with an emission current of 
100 yamps, beam current of 0.5 yamps , and an operating 
voltage of 25 kv. The intensities of the gold and silver 
La lines were measured with a LiF and ADP crystal, respec- 
tively. Two (2) measurements for a fixed time count of 10 
seconds were taken at five different locations on each film. 
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Background measurements were taken by moving the spectrome- 
ters off the peak position. 

Thin - Film Couple (Ag- Au) 

Oxidized Silicon wafers with the Ta annulus sputtered 
onto the surface were used as substrates. The substrates 
were mounted onto the substrate holder and clamped firmly 
into place by the copper restraining plates. A Personna 
74 double edge razor blade, with its commercial stamping 
removed, was cleaned by ultrasonic agitation in reagent 
grade acetone and trichloroethylene. A small strip of 
.009 mm aluminum foil was spot welded onto the razor blade 
to cover the center slot design. The final assembly was 
cleaned further in solvents and then placed into the razor 
blade support mechanism. The substrate holder was placed 
into the vacuum system. The source barrier was then 
adjusted horizontally to be true to the razor blade edge, 
and vertically so that as the substrate was heated and 
expansion occurred, the vertical height between the top of 
the barrier and the razor blade edge was as small as 
possible . 

The vacuum system was pumped down to 5 x 10" 7 torr while 
the substrate was heated by the variac to the annealing 
temperature. The sources were preheated as described above. 

While the substrate was held above room temperature, 
helium quench gas was slowly passed through the quench tubes. 
The small flow rate was monitored by terminating the gas 
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line in a flask filled with distilled water. This prevented 
oxidation of the inner surfaces of the copper quench tubes, 
which were exposed to 1 atm. pressure. 

When the substrate was at the annealing temperature, 
the Leeds and Northrup proportional counter was switched 
into the heating circuit, controlling the temperature of 
the substrate to within + 1°C. 

Deposition of the gold and silver, while the substrate 
was held at the annealing temperature, was carried out simul- 

O 

taneously. The rates (5 A/sec) of each source were adjusted 
equally by the crystal monitors, prior to opening the 
shutter, so that the time of deposition and film thickness 
of the gold and silver were the same. Silver and gold of 
99.99% and 99.9999% purity, respectively, were used as 
source material for deposition. 

At the conclusion of the film deposition, the anneal 
was begun. The anneal was terminated by removing the flask 
from the terminus of the gas line and then forcing helium 
gas at a pressure of 12 lbs./sq. inch through the quench 
tubes. The cooling curves were exponential in character, 
representing quench rates of approximately 75-100°C per 
minute, down to substrate temperatures of around 60°C. The 
quench rate lowered to approximately 50°C per minute for the 
two low temperature anneals. 

The substrates were then removed from the vacuum system 



for further analysis. 
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Film thickness of both the pure gold and pure silver 
portion of the couple and thickness variation across the 
region of mixing were measured by Dektac. 

X-ray emission data were taken with the ARL-AMX 
Electron Microprobe. Intensity measurements were obtained 
with the electron microprobe operating as described above 
in the subsection on thin film standards. Intensity profiles 
as a function of distance were determined by recording the 
characteristic x-ray intensity while the sample was moved 
under the electron beam at a rate of 4 microns/minute. 

Fixed time counts of 10 seconds were taken as the sample 
was scanned. 

Photomicrographs of the film surface were taken on a 
Leitz orthoplan microscope and a JSM-U3 Scanning Electron 
Microscope. The scanning electron microscope was operated 
at 27 kv with a sample current of 1.89 x 10" 10 amps. 

Limited x-ray emission data were also taken with the scanning 
electron microscope, taking advantage of the small spot size 
of the electron beam, and will be noted as the data are 
reviewed . 

Grain sizes of the films were determined using the 
transmission mode of a Siemens Elmiskop 101 electron 
microscope operating at 100 kv, with beam currents of 10 
and 60 microamps. 
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Incremental Couple (Au- 201 Ag ) 

Oxidized Silicon wafers with the Ta annulus sputtered 
onto the surface were again used as substrates. The thin 
film deposition was carried out in two steps. 

Initially the substrate, mounted onto the substrate 
holder, was placed into the vacuum system without the razor 
blade. The substrate was heated to the annealing tempera- 
ture, and then a layer of gold film was deposited across 
the entire substrate surface. After cooling to room tempera- 
ture, the substrate holder was removed from the vacuum sys- 
tem where the razor blade was mounted into the support 
mechanism. After alignment of the razor blade with the 
barrier mechanism, as in the case of the Ag-Au couple, the 
system was again evacuated. The substrate was again heated 
to the annealing temperature, and films of gold and silver 
were deposited over the initial gold underlayer, in the 
manner of the Ag-Au couple described in the previous sub- 
section. 

Film thickness, x-ray emission data, photomicrographs 
and grain size determinations were taken in the same manner 
as described in the previous subsection. 
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RESULTS 

Thin Film Standards 

The determination of the thin film empirical conversion 
parameter between the measured relative intensity of the 
microprobe and thin film composition, involved the fabrica- 
tion of twenty-four (24) different films. Gold-silver films 
of 0, 20, 40, 60, 80, 100 atomic % silver were each made of 
4 different approximate thicknesses, 500, 1,000, 1,500 and 
2,000 X. Each film thickness was measured on the Dektac 
where the accuracy of measurement is ± 50 X. Two plots of 
measured characteristic x-ray intensity versus measured 
film thickness, were made for each film composition. The 
two plots were of the characteristic x-ray data of the La 
lines taken from the gold and silver spectrometers. A 
linear least squares fit was made on each intensity versus 
thickness plot. 

From the resultant linear analytical equation for each 
plot, the film x-ray intensity was determined for numerical 
thicknesses of 500, 1,000, 1,500 and 2,000 A. K^g and K^ u 
for each composition and film thickness, in the form of 
equation (1) , were determined after correcting both the 
pure and alloy films for their respective backgrounds. 
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The backgrounds of the pure Au and Ag films were deter- 
mined by moving the spectrometers off each side of the 
characteristic La peak of the Au and Ag line respectively, 
and averaging the two values. The alloy background was 
determined by the following relationship: 



B 1 1 = B . C. 

alloy pure Au Au 



+ B . C . 

pure Ag Ag 



(15) 



where B a n 0 y = t ^ ie background of an alloy film of 

thickness, t. 

Bpure a = t ^ ie background of a pure A film of 
thickness, t. 

C = atomic fraction. 

The values of , K^ u and , C^ u respectively, were 
plotted in the form of equation (3) , where another linear 
least squares fit was made, the slope representing the value 
of the empirical conversion parameter, & A g Au and a AuA g 
respectively. The process was repeated for each film 
thickness. In addition, a linear least squares fit was 
made for all the data of each spectrometer to determine 
an average value of the empirical conversion parameter for 

O 

film thicknesses between 500 and 2,000 A. 

The statistical calculations discussed above were 
accomplished by the computer program in Appendix C. Tables 
1 and 2 are the direct computer output printed on the Gould 
electrostatic printer, displaying the results of the cali- 
bration data taken off the silver and gold spectrometer 
respectively . 
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] ncremontal Couple ( 80%Au- Au ) 

Table 3 summarizes the experimental runs completed 
during the study of chemical diffusion in the incremental 
couples. The couples consisted of a pure gold film layered 
side by side to an alloy film of equal thickness, and with 
a composition of 80 and 20 atomic % gold and silver, 
respectively. The film thickness for all the runs was 

O 

approximately 1,500 A; actual measurements made on each 
film using the Dektac are reported in Table 3. 

Five (5) couples were annealed for 24 hours at tempera- 
tures of 200, 250, 300, 350 and 416°C. Two additional 
couples were annealed at 416°C for 12 and 48 hours. 

One (1) couple (Run 227) was not annealed and represents 
the initial (time = 0) concentration profile. Figure 16 
shows the concentration-distance profile generated from the 
computer plotting program of Appendix D. The computer pro- 
gram reads in the x-ray intensity data from the microprobe 
as a function of lateral distance (Stepsize 0.8 microns), 
calculates the measured relative intensity corrected for 
background, and converts K^g to units of mass concentration 
(grams Ag/cc) through the empirical conversion parameter 
and appropriate value of alloy density. 

Figures 17 through 23 show the concentration-distance 
profiles, generated from the same computer program, of the 
annealed couples. In order to compare the annealed profiles 
with the initial profile, and thus determine values of 
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Table 3. Summary of Experimental Conditions for Runs 224-232 





Incremental 


Couples (80Au-Au) 


Run 
No . 


Temperature 
of Anneal 
(°C)* 


Time 

of Anneal 
(hrs) 


Film , 

Thickness 7 Comments 

(A) 


224 
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1500^ No interface contact 
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300 


24.0 


1480 
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0.0 
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1480 


229 


416 


48.0 


1440 


230 


250 


24.0 


1525 


231 
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12.0 


1530 


232 


200 


24.0 


1500 


* + 


1°C 






t + 


50 A 






1 not measured 
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chemical diffusivity as a function of temperature, the 
following analytical procedure was used. 

In the case of the thin film incremental couples, it 
was assumed that Pick's second law for one dimensional 
diffusion is a valid approximation to the measurement of the 
rate of change of concentration with respect to time, and 
is written; 



3C _ n a 2 c 
3t " U 8 x 2 



(16) 



where D f D(C) ; D f D(x) 

C = mass concentration 
t = time 
x = distance 

The solution to equation (16), if we have a thin 
respect to the lateral direction) strip centered at x 
with an initial thickness (with respect to the lateral 
direction) 6 and concentration Cj is 



C(x , t) = Cl ( 0 - ? 0),S exp 
2/irDt 



x‘ 

"3Dt 



(wi th 
= 0 , 



(17) 



A thin strip centered at x = 6, with an initial thick- 
ness 6 and concentration C 2 , gives 



C (x , t) 



C 2 (6,0)6 

2/iUt 



exp - 



(x-6) 2 ] 
TDt j 



(18) 



For an arbitrary collection of thin strips of thickness 6, 
for example, and arbitrary initial profile 
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c ( x >t) = I t' n (x,t) 



+ oo 




The thin film incremental couple was assumed to 
represent a collection of thin strips of lateral thickness 
6, and perpendicular thickness t (measured film thickness). 
The lateral strip thickness 6 was chosen to be 0.8 microns, 
the distance between experimental data points. For the case 
of the initial concentration profile of Figure 16, the con- 
centration of each thin strip was determined by the follow- 
ing relationship 



or the average concentration between two successive experi- 
mental concentration data points. 

A computer program was written to numerically solve 
equation (19), for time = 24 hours, and a series of 
expected values of D, using the experimental data points 
of the initial concentration profile, and is shown in 
Appendix E. 

On the basis of the numerical results of the computer 
solution to equation (19) the linear least squares slope of 
the numerically solved concentration profile, between the 
values of 1.56 and 0.56 grams Ag/cc, was chosen as the test 
parameter to compare the experimental profiles of Figures 
17 through 23. 



C n (n6,0) = C n (n6,0) + C n+ x ( (n+1) 6 , 0) /2 



( 20 ) 
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The linear least squares slope was determined using the 
statistical procedures of Volk [39], and included a test 
for the standard error of the estimate, sum squares deviation 
of y (concentration), coefficient of correction, and variance 
of slope. The statistical procedures are written into both 
the computer plotting program for the experimental profiles, 
(Appendix D) , and the computer program which solves equation 
(19) (Appendix E) . 

Table 4 summarizes the statistical calculations for the 
annealed profiles and compares them with the numerical solu- 
tions of equation (19) . 

From the comparison of slopes of the annealed profiles 
with the numerical slopes of the solution of equation (19), 
values of D were determined. The concentration profiles 
generated from the numerical solution of Run 227 using equa- 
tion (19) are drawn for the appropriate value of D on each 
experimental concentration profile shown in Figures 18-23. 
Figure 24 shows a plot of In D versus 1/T for the 24 hour 
anneals at 200, 250, 300, 350 and 416°C. The experimental 
data fit the equation 



D = 1.9x10 



exp 



10,900 

RT 



cm 2 sec - 1 



( 21 ) 



In addition, the values of D, for the 12 and 48 hour 
annealed profiles at 416°C, are reported in Table 4. These 
values were determined in a manner similar to the 24 hour 
anneal, with the slope of their respective concentration 
profiles between 1.56 and 0.56 grams Ag/cc compared with 



Table 4. Summary of the Statistical Calculations of the Incremental Thin Film Couples. 
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°C 




Figure 24 

Plot of Ln D versus 1/T for the Incremental Thin Film 
Couples Annealed for 24 Hours. 
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the numerical solution with the corresponding times entered 
into equation (19) . 

Grain size determinations, using the transmission mode 
of the electron microscope, were made on each couple. The 
grain size varied between 1,000 and 2,000 A on all couples 
observed. No variation across each couple, including the 
diffusion zone was observed. Figure 25 shows photomicro- 
graphs of the grain distribution in the 200 and 350°C 
couples. Figures 26 through 28 are photographs of the 
surfaces of selected incremental couples showing the change 
in surface structure with increasing deposition temperature. 

Ag-Au Couples 

Table 5 is a summary of the experimental test conditions 
of Runs 201 - 223. Each run involved the deposition of a 
pure gold film layered side by side with a silver film of 
equal thickness. 

Runs 202 through 212 have comments in Table 5, listing 
reasons for the author's inability to use them for quantita- 
tive results. The problem occurring in the initial runs 
(202 - 206) was film breakup (dewetting) of the gold side of 
the couple. Along with this, there was some interesting 
particle shaping of the resultant discrete gold islands 
observed, such as smooth equilateral triangles, like those 
observed [40] in epitaxial gold films on [111] surfaces of 
M 0 S 2 , and particle shapes observed by Sundquist [41]. 
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Figure 25 

Transmission Electron Micrograph of Section of the Surface 
of Unsupported Thin Film. 23, 000* 

A. Gold Side of Incremental Thin Film Couple 

Annealed at 200°C for 24 Hours. 

B. 80%Au-20%Ag Side of Incremental Thin Film 
Couple Annealed at 350°C for 24 Hours. 
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SOfcu - 2oft 9 




Au 



Figure 26 

Scanning Electron Photograph of the Surface of the 
Unannealed Incremental Thin Film. Couple (Run 227) Showing 
Interface Region Film at 45° Tilt to the Electron Beam 

10 , 000 *. 
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Figure 27 

Scanning Electron Photograph of the Surface of the 
Incremental Thin Film Couple Annealed at 250°C for 24 Hours, 
Showing Interface Region. Film at 45° Tilt to the 



Electron Beam 



10 , 000 * 
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Figure 28 

Scanning Electron Photograph of the Surface of the 
Incremental Thin Film Couple Annealed at 416°C for 24 Hours, 
Showing Interface Region. Film at 45° Tilt to the 
Electron Beam. 10, 000*. 
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l’;i b 1 c 5. Summary of experimental Conditions of Runs 201 
through 223. (Ag-Au) Couples 





Temperature 
o f 


Time 

of 


Film 






Anneal* 


Anneal 


Thickness 




Run 


(°C) 


(hrs) 


(A) + 


Comments 


201 


25 


0.0 


1000+ 


used as time = 0 profile 
for Ag-Au couples. 


202 


675 


5.0 


1000+ 


excessive breakup 
(dewetting) of gold 
portion of couple 


203 


575 


16.00 


iooot 


same as 202 


204 


500 


2.0 


1000 + 


partial dewetting of 
gold side; no interface 










contact 


205 


675 


0.0 


1080 


excessive surface rough- 
ness of gold and silver; 
partial dewetting 


206 


500 


0.0 


1000 + 


same as 205 


207 


428 


0.0 


1000 + 


no interface contact 


208 


416 


24.0 


1550 


wide region of mixing 


209 


416 


48.0 


1500 + 


no interface contact 


210 


416 


48.0 


1500 + 


same as 209 


211 


416 


48.0 


1500 + 


same as 209 


212 


416 


48.0 


1500 + 


same as 209 


213 


416 


48.0 


1000 


run 201 annealed 


214 


416 


96.0 


1000 


run 201 annealed; 
dewetting of gold 


215 


416 


72.0 


1000 


run 201 annealed 


216 


416 


24.0 


1520 




217 


416 


48.0 


1540 




218 


416 


12.0 


1510 




219 


416 


2.0 


1490 




220 


300 


48.0 


1560 




221 


300 


24.0 


1500 + 


no interface contact 


222 


300 


24.0 


1470 




223 


300 


12.0 


1490 





* t 1°C 0 
t t 50 A 

•f not measured on Dektac 
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Three (3) runs, 213, 214 and 215 were made, using 

O 

sections of the 1,000 A unannealed couple (Run 201), at 
416°C for times of 48, 72 and 96 hours. The formation of 
a porous zone, running parallel to, and on the silver side 
of the interface, was observed in the couples annealed for 
48 and 72 hours. In addition, a zone of protrusions was 
observed on the gold side of the interface. The porous 
zone was so extensive that complete film separation occurred. 
The width of the porous zone was slightly larger on the 72 
hour couple. In the couple annealed for 96 hours, the film 
separation due to the porous zone was evident. However, 
dewetting was observed on the gold side of the film, with 
particle shaping like that described above for films 
deposited at higher temperature. The silver side of the 
couple also exhibited extensive pitting, and investigation 
on the scanning electron microscope revealed a high degree 
of surface roughness. Subsequently, no quantitative study 
of diffusion was undertaken for these runs. 

O 

Runs 216 through 219 involved the deposition of 1,500 A 
thin film couples at 416°C, for annealing times of 2, 12, 

24 and 48 hours. Again, the formation of the porous zone 
and surface protrusions were observed on the silver and 
gold side of the original interface, respectively. Figures 
29 through 32 show photomicrographs of the respective inter- 
face region of the four (4) couples. Figure 33 is a scanning 
electron photograph of the interface region showing the 
porous zone, and the associated surface protrusions. 
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Figure 29 

Optical Microscope Photograph of Interface Region of Ag-Au 
Thin Film Couple Annealed at 416°C for 2 Hours. 1080x. 
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Figure 30 

Optical Microscope Photograph of Interface Region of Ag-Au 
Thin Film Couple Annealed at 416°C for 12 Hours. 1080x. 
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Figure 31 

Optical Microscope Photograph of Interface Region of Ag-Au 
Thin Film Couple Annealed at 416°C for 24 Hours. 1080x. 
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Optical Microscope 
Thin Film Couple 



Figure 32 

Photograph of Interface Region of Ag-Au 
Annealed at 416°C for 48 Hours. 1080x. 
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Figure 33 

Scanning Electron Photograph of Interface Region of Ag-Au 
Thin Film Couple Annealed at 416°C for 24 Hours. 

Film at 45° Tilt to the Electron Beam. 3,000x. 
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An investigation of the chemical composition of the 
surface protrusions was undertaken by measuring characteris- 
tic x-ray lines generated from the electron beam of the 
scanning electron microprobe. For a counting time of 200 
seconds, the counts of the La lines of the silver and gold 
remained constant with the electron beam, on and off, 4 
different surface protrusions. The silicon counts decreased 
to one-half their original value when the electron beam was 
on a surface protrusion. 

Characteristic x-ray measurements were also taken for 
Runs 216 - 219 as a function of lateral distance on the 
electron microprobe. A computer program similar to that 
used for the incremental couples, converted the intensity 
counts to units of mass concentration. The direct computer 
output of the couples are shown in Figures 34-37. 

Because of the formation of the porous zone and associ- 
ated protrusions, there was film thickness variation as a 
function of lateral distance, thus the concentration-distance 
profiles were suspect. 

A scanning electron microprobe investigation of the 
couples, however, indicated the presence of silver on the 
gold side of the porous zone extending away from the inter- 
face approximately the same distance as indicated by the 
concentration-distance profiles. No counts of gold were 
observed on the silver side of the porous zone. 
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The dip in the concentration value observed in the 
concentration-distance profiles was taken as a representa- 
tive of the center of the porous zone. A computation of the 
area under the concentration profile from the center of the 
dip on into the gold region (increasing relative distance) 
was made on each profile. The area was determined using a 
computer program written to numerically integrate, employing 
Simpson's rule. The results are summarized in Table 6. 

A computer program similar to that of Appendix E was 
written to treat the peak in concentration observed in the 
4 couples. Run 219 (Figure 34), the couple annealed for 2 
hours, was treated as the initial profile (time = 0) and 
the values of concentration on the gold side of the porous 
zone was fed into the computer program along with a compli- 
mentary set of image concentration values to make the initial 
profile an even function. The reduction in peak height of 
the 12, 24 and 48 hour annealed couples was compared with 
the numerical solution for times of 10, 22 and 46 hours 
respectively. The numerical solution to the concentration 
"pulse" of Run 219 is drawn on each concentration profile 
for an appropriate value of D. D values of 1 , 5 and 9xl0 -11 cm 2 / 
sec fit the experimental profiles for the 12, 24 and 48 hour 
anneal, respectively (Figures 35-37). 

Runs 220 - 223 were 1,500 A thick thin film couples of 
Ag and Au annealed at 300°C for 12, 24 and 48 hours. The 
formation of a porous zone and associated surface protrusions 
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able 6. 



Run 

219 

219 

218 

218 

216 

216 

217 

217 



Comparison of Area under Concentration Curves from 
the Edge of Porous Zone into the Pure Gold Region 
of the Couple. For Runs 216-219, Annealed at 416°C 



Time of 

Microprobe Anneal 
Scan No. (hours) 



Area 
(cone . 
microns) 



Average 

Area 



1 2 61.9 58.4 

2 2 54.8 



1 12 62.1 60.2 

2 12 58.4 



1 


24 


54.0 


2 


24 


49.3 



1 


48 


40.1 


2 


48 


35.3 



51.7 



37.7 
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were again evident. The width of the porous zone of the 12 
hour annealed couple was similar to that of the 2 hour 
annealed couple at 416°C. The porous zone enlarged with 
time, and film separation was evident in the couple annealed 
for 48 hours (Run 220). X-ray measurements as a function 
of lateral distance were made on the electron microprobe 
and the data were converted to units of mass concentration 
as described above. However, the porous zone of the 12 hour 
annealed sample was so small in extent that the dip in con- 
centration was not evident. Therefore, the semi-quantita- 
tive technique to determine the rate of silver penetration, 
as described above, was not attempted. 
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DISCUSSION 



Thin Film Standards 

The linearity (high degree of correlation) observed 
between the microprobe intensity and film thickness for all 
the data exhibited in Tables 1 and 2 suggests that the thin 
film approximation, intensity is proportional to mass thick- 
ness, is a valid assumption throughout the entire composition 
range of the gold-silver system, for the film thickness 
studied. The reader is cautioned, however, that the numeri- 
cal measurements are only valid for gold and silver films on 
a silicon substrate with a 1,500 A oxidized layer on the 
surface . 

The statistical equations used in the analysis of thin 
film standards, and subsequent work on thin film couples, 
are shown in Appendix F. 

The maximum deviation from linearity occurred with the 



2,000 A films as should be expected, especially in the case 
of pure gold films where I/I^ulk was approximately 0.6. The 
maximum relative error, .25, in the empirical conversion 
parameter, was observed in the measurements taken off the 



silver spectrometer for films between 500 and 2,000 A. 
Because of the small relative error, one conversion param- 
eter was used. The parameter was determined from a linear 
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least squares fit of all the data, representing the relation- 
ship between measured relative intensity and weight fraction 
composition for films of thickness between 500 and 2,000 A. 
Again, the largest deviation from the straight line fit 
occurred for the 2,000 A films of high gold composition, 
where the thin film approximation became suspect. 

The relation between an error in the empirical conver- 
sion parameter and errors in the calculated composition value 
is given by the usual expression 

AC = g da (22) 

and from equation (2) 

= C(l-C)/ a . (23) 

The error in composition is 

AC = C (1 - C) Aa/a (24) 

The maximum absolute error in composition is one-fourth (£) 

of the relative error in the conversion parameter, or 
(.25) (.25) which is approximately 6 weight %. 

This, however, is based on the assumption that the com- 
position of the standards is known to within negligible 
limits of the error of the calculated value. 

The determination of the composition of the thin film 
standards used in the experimental work of this author was 
based on the assumption that the thin film density in both the 
pure and alloy states was that of the reported bulk value 
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discussed in the Background and Literature Survey section. 

Hartmann [42] obtained bulk values for Au and Ag at 
thicknesses down to 200 A. However, Blois and Reiser [43] 
found the density of Ag films to approach bulk values at a 
thickness of 1,200 A. The data of [43] are suspect [57], 
for silver, which readily forms a surface oxide layer, 
lowering the effective density of the film. 

This author, in the initial stages of determining the 
orientation relationship between the mass deposited on the 
monitoring crystals and that of the substrate, made three 
determinations each between the frequency shift of the 
crystal and measured thickness (Dektac) of pure gold and 
silver films. Using the relationship, 

Af = Bpt > (25) 

where Af = frequency shift of the crystal 
t = film thickness 

0 = substrate-crystal-source orientation parameter, 
p = film density 

a linear relationship was observed between Af and the 
measured value of film thickness. Solving equation (25) 
for and p^ u respectively, a bulk value was obtained 

within the relative error of the measured film thickness. 

If the standards are assumed to be of bulk density, then the 
major source of error is the measurement of film thickness, 

+ 50 X. The error in composition can be written 
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At. / t A 

AC (atomic %) = ^ — S — — (26) 

Ag + Au 

t Ag t Au 

which has a maximum value of ± 10 atomic % for film thickness 

O 

of 500 A and at each extreme of the composition range. 

Ag- Au Couples 

The experimental observations of film breakup on the 
gold side of the thin film couple suggest that the mechanism 
of atomic reorientation manifests itself as a function of 
time and temperature, and can be attributed to an annealing 
effect. The agglomeration of the gold side of the couple 
observed in Runs 202 - 206, at substrate temperatures above 
500°C, was only observed at the conclusion of the anneal. 

The initial condition of the film, immediately after film 
deposition was in doubt and the observed effects, when com- 
pared to that of Run 201, could only be attributed to the 
substrate temperature. However, in the series of runs at 
416°C (Runs 213 - 215) the 1,000 A film deposited at room 
temperature (Run 201) was used, and the following anneal led 
to agglomeration of the gold after 96 hours (Run 214) . 

If a polycrystalline material is annealed, grooves 
develop at the intersection of the grain boundaries with the 
metal-vapor surface. If it is assumed that the grain bound- 
aries do not move during annealing, then it can be expected 
that the groove depth (d) deepens with time. Mullins' [44] 
well known analysis of grain boundary grooving kinetics was 
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applied to the thin films of the present work, thus 

(0.486)y b 



d = 



D y c a 2 NL 

s 1 S t 

FT 



(27) 



where d = depth of groove 

Yg = specific free energy of the grain boundary 
Yg = specific free energy of the surface 
D s = surface diffusion coefficient 
k = Boltzmann's constant 
T = temperature of anneal (°K) 
t = time of anneal 

N = number of surface atoms/unit area 
Q = volume per atom 

In order to relate the experimental observations of the 
gold film breakup, a two stage model was proposed. In the 
initial stage there would be groove deepening to the sub- 
strate, the depth of the groove equal to the original film 
thickness. Secondly, the hole grows, with critical size 
agglomerates (adjacent to the holes) growing at the expense 
of the surrounding metal. 

If equation (27) is evaluated at 416°C, (temperature of 
anneal of couple where gold agglomeration occurred after 96 
hours) using the following values, 

from [49] 
from [49] 



Y s = 1450 



Y B = 450 

12 = 1.7 x 10" 2 3 cm 3 



N = 1.2 x 1.15 cm' 
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D g = 5 x 10~ 8 cm 2 /sec for Au at 416°C from [51] 

1> S = 6 x 10" 11 cm 2 /sec for Ag at 416°C from [51] 

The depth of the groove would equal 1,000 A (film thickness) 
for annealing times of 62 and 2.5 * 10“* hours for gold and 
silver respectively. 

In the case of gold, we would expect groove deepening 
to the substrate after 62 hours. Following this, the anneal 
enables the film substrate system to lower its free energy 
by lowering its interfacial free energy. The motive force 
for the growth of an existing hole by surface diffusion is 
the surface tension vector sum of the film-substrate, sub- 
strate-vapor and film-vapor interfaces. The flow of migrating 
atoms is due to the difference in chemical potential at regions 
of differing curvature. 

Presland et al . , [45 and 46] in a series of two articles 
applied the theories of Mullins [44] to his experimental 
observations in annealed thin films of silver. He showed 
that following an induction period (hole deepening to the 
substrate with subsequent hillock growth) that further hillock 
growth followed, with critical size agglomerates growing at 
the expense of surrounding areas. The continuing hole growth 
resulted in a decrease in metal-substrate interfacial area, 
resulting in the separation of agglomerates (film breakup) . 

Brandon and Bradshaw [47] showed that the surface self 
diffusion coefficient of silver is increased in the presence 
of oxygen. Presland could only produce agglomeration in his 
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high vacuum deposited silver films when the silver vapor- 
metal interface was exposed to oxygen pressures greater than 
I0' 3 torr. Below p n = 10” 3 torr, he reported groove deepen- 
ing with subsequent hillock growth. However, the catchment 
area for the critical size agglomerates increased rapidly 
with decreasing Pq 2 > suggesting changes in surface tension 
and the surface self diffusion coefficient with oxygen 
partial pressure. In comparison, gold agglomeration has been 
observed [48] in unsupported thin films in high vacuum. 

The equilibrium shape of a particle on a substrate is 
determined by the vector sum of the interfacial free energies 
of the film, vapor, substrate system. If the observations of 
[46 and 47] are assumed to be operative for the gold film, 
then in conjunction with the sufficient groove deepening, as 
predicted by equation (27) , the vector sum of the interfacial 
free energies is in accord to predict hole growth and subse- 
quent agglomerates. 

In the case of silver, the reduction of the film-vapor 
interfacial free energy, due to the presence of oxygen on 
the free surface of the silver as suggested by Presland [45, 
46], would not be sufficient cause for hole growth. 

More likely, the increase in the surface self diffusion 
coefficient due to oxygen pressure would reduce the time 
necessary for groove deepening to the substrate. Following 
groove deepening, a reduction in the substrate-vapor inter- 
facial free energy due to oxygen pressure would be favorable 
for hole growth followed by film agglomeration. 
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The observed faceting on the gold film agglomerates 
could be attributed to the reported [41] surface anisotropy 
of gold. Gjostein [50], who provided some of the impetus for 
the studies of [45 and 46] in an earlier paper [51], is 
currently developing a manuscript to possibly improve on the 
studies of [45, 46 and 47] and should provide interesting 
reading when released. 

To avoid the agglomeration of the gold film, further 
runs (216-223) were made at temperatures below 500°C and for 
annealing times only up to 48 hours. In each of these runs 
the porosity observations of bulk investigators with pure 
silver and gold sandwich couples were similarly reproduced 
in the thin film couples. The porosity increased with time 
and temperature of anneal. However, in contrast to the bulk 
counterpart, physical separation of the films occurred in 
the lateral direction. 

The physical separation of the two regions of the film 
leads to thickness variations of the couple around the porous 
zone. The concentration profiles of Figures 34 to 37 were 
determined by assuming a uniform thickness across the diffu- 
sion zone. Therefore, in the region of the porous zone, the 
concentration profile is in error. 

Secondly, the physical separation of the couple severely 
limits the atomic mobility of the two species. When this 
occurs, the quantity of gold and silver which had previously 
diffused into the silver and gold respectively, remains fixed 
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with time. The observation of the similarity of the areas 
under the concentration profiles of the 2, 12 , 24 and 48 hour 
anneals confirms this conclusion. The smaller area observed 
in the 48 hour anneal (Run 217) can be attributed to the 
inability of the electron microprobe to detect the small 
concentration of Ag above background at distances far from 
the porous zone, and the corresponding increase in error of 
the empirical conversion parameter with low values of solute 
concentration . 

The mass of silver remaining fixed with respect to time, 
after the formation of the porous zone, creates a set of 
boundary conditions different from that of an initial t = 0 
profile. The series of photographs of the four anneals shown 
in Figures 29-32 indicate that physical separation of the 
couple is almost complete after the 2 hour anneal (Figure 29) . 
On either side of the porous zone the film thickness is 
uniform, therefore the concentration profile of the 2 hour 
anneal can be taken to represent the new initial profile 
(t' = 0) and the 12, 24 and 48 hour anneals to represent 
values of t' = 10, 22 and 46 hours respectively. 

The thin strip approximation was applied to the concen- 
tration "pulse" of silver observed on the gold side of 
porous zone of the t' = 0 profile. The center of the "pulse 
peak" indicated on Figure 34 was taken to represent the edge 
of the porous zone shown in Figure 29. 
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The exact solution [52] to the initial (t' = 0) concen- 
tration profile would be to represent the profile by f(x) 
and solve using a heat flow analogy. The porous zone would 
represent an impenetrable barrier (insulating medium) . 

A numerical approximation to the exact solution was used, 
incorporating the thin strip approximation. The concentra- 
tion "pulse" was made into an even function by reading into 
a computer program the image concentration values as shown 
in Figure 34. The thin strip approximation predicts a 
reduction in peak height with points of constant concentra- 
tion moving laterally parabolically with time. The resultant 
fit to the t' = 10, 22 and 46 hour profiles shown in Figures 
35 to 37 lends validity to the above approximation. 

The reduction in D as the time of anneal increases can 
again be attributed to the errors in detecting small solute 
(Ag) concentration values and the rising error in the con- 
version parameter with decreasing solute concentration. 

Also, the errors in the thin strip approximation decrease 
as /Dt increases relative to the assumed (0.8 microns) strip 
thickness (lateral direction). The value of D therefore, 
is probably somewhere between the results for the t' = 10 
and 22 hour anneals. Averaging the two results gives a 
D of 7 x 10 -11 cm 2 /sec. which closely represents what was 
reported for the incremental couples at 416°C. 

An interesting observation was that of the secondary 
zone of porosity observed in the photographs for the 24 and 
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48 hour anneals (Figures 31 and 32) . This could possibly be 
attributed to the new set of boundary conditions operative 
after the 2 hour anneal. The mass flow of silver in the t' 
framework could be viewed as a thin film of gold layered 
side by side of an alloy film of say 50 atomic % silver with 
only a 30 micron extent (f width of concentration "pulse") 
in the lateral direction. The new set of boundary conditions 
would also shift the position of the original interface and 
the corresponding position of the maximum vacancy flux 
gradient [53] . 

This reasoning implies that the vacancy flux gradient 
must be a strong driving force for the nucleation of voids 
by an assimulation of excess vacancies. On the basis of 
energy considerations it would appear logical for the 
vacancies accumulating in the subsequent mass depletion zone 
to continue to annihilate themselves on the already formed 
voids in the original porous zone. 

The x-ray investigation on the scanning electron micro- 
scope of the surface morphology on the gold side of the diffu- 
sion zone confirmed the suspicion that the morphology repre- 
sented positions of increased film thickness or surface pro- 
trusions . 

As with the porosity formation, no in-depth investigation 
of the formation mechanism of these surface protrusions, or 
if the reader will allow, surface hillocks, was undertaken, 
for it was not the purpose of this work. 
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However, some speculations will be made from the previ- 
ous work of bulk investigators in the Ag-Au system and those 
studying thin film annealing phenomenon. 

Several investigators, [54, 55 and 56], studying thin 
film annealing phenomenon have reported hillock growth in 
thin films, including Au , and attributed it to a mechanism 
of stress reduction. The driving force in their studies 
was the differential stress between the film and substrate, 
developed by the AT between the deposition and annealing 
temperature. In the present work, deposition at the anneal- 
ing temperature (excepting Runs 212 - 215) eliminated this 
as a driving force. However, Barnes [30] and Balluffi [32] 
reported strain in the diffusion zone of bulk sandwich 
couples of the Cu-Ni and Ag-Au systems respectively. The 
reduction of strain in the diffusion zone and/or dislocation 
climb creating planes of vacancies similar to that proposed 
by Ruth [35] in the silver-gold system, could account for 
the observed hillock growth. 

Incremental Couples ( 80 Au-Au ) 

As in the bulk studies of Johnson [38] with sandwich 
couples with small initial concentration differences, no 
porosity or hillock formation was observed in the diffusion 
zone of thin film couples with an initial 20 atomic % con- 
centration difference. 
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The transmission electron microscope studies of the 
couples' grain size led to the conclusion that the grain size 

O 

(1,000 - 2,000 A) of the films was independent of the anneal- 
ing temperature, and film composition over the working range 
of the couples. The grain size of Ag films deposited at 

O 

room temperature has been reported to be as low as 50 A. 
However, with the high deposition temperature of the present 
study, the grain size is expected to be larger, due to the 
increasing surface mobility of adatoms and clusters during 
deposition . 

The surface roughness, which was uniform over each 
individual couple, increased with annealing temperature as 
was shown in Figures 26 to 28. As the surface mobility of 
adatoms increases (higher substrate temperatures) , the con- 
densation can occur preferentially at the surface concavites , 
and thus tend to smooth the surface. However, as was 
reported [41] with gold and experimentally observed in the 
present study, anisotropy of the surface energy will cause 
the condensate to grow preferentially along certain crystal 
faces, enhancing the resultant surface roughness. This is 
strongly evident in the couple annealed at 416°C shown in 
Figure 28. 

In each of the concentration profiles, the variance in 
the slope was within .005. The value of .005 was within the 
statistical error of each slope. The error was determined 
from the experimentally observed variation in slope between 
successive scans of the microprobe across the diffusion zone 
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on different portions of the couple. The three values of D 
determined at 416°C for annealing times of 12, 24 and 48 
hours lent credence to the use of the thin film approximation 
which predicts parabolic movement of atomic planes of 
constant composition with time. 

The values of D determined for the 24 hour anneals at 
temperatures of 200, 250, 300, 350 and 416°C exhibited 
linearity on a In D versus 1/T plot with correlation 0.99. 
This indicated a single atomic transport mechanism to be 
dominant over the working temperature range of the couples. 
The value of activation energy was almost £ of that reported 
for bulk diffusion in Ag-Au couples. 

In the thin film couple, three possible atomic transport 
mechanisms are considered to be operative: Transport of 

matter along the grain boundaries and the substrate-film 
interface, through the crystal lattice, and along the film- 
vapor interface. The apparent diffusivity can be written as 



U 



app 



D GB 8 + D L d-g-s) * D s s 



(28) 



where 



D 



GB 

D, 



D 



s 



grain boundary diffusion coefficient 
lattice diffusion coefficient 
surface diffusion coefficient 

fraction of time a diffusing atom spends on a 
grain boundary or along the substrate- film 
interface 

fraction of time a diffusing atom spends on the 



surface . 
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For the case of a binary alloy, all atoms are not equiv- 
alent so the total time a diffusing atom spends on a grain 
boundary or surface is equal to fractional number of sites 
which are in or on each respectively, times its atomic 
fraction . 

To estimate g, consider a sample made up of "spheres" 

that touch each other completely. The total number of atoms 

in a "sphere" is 4/3irr 3 /Q, where Q equals the atomic volume. 

The total number of atoms in the grain boundary equals 
2/3 

4frr 2 /ft 7 . Therefore, 



g 



(atm. 



frac . ) 



' 4ur 2 


to 


^ 2/s J 


4nr 3 

V / 



(atm. 



frac . ) 




l/3 0 

for 7 = 2.2 A 

r = 1,000 A 
atm. frac. = 0.1 



Therefore g = 10 3 

To estimate s, consider a surface monolayer of thickness 

O O 

2 A. For a film thickness of 1,500 A and atomic fraction 
0.1, s = 10' 4 . 

We can rewrite equation (28) as 

D app ' 10 ' 3 D GB * D L * 10 "' D s < 29 > 

The activation energy for chemical diffusion through the 
lattice is approximately 40 kcal/mole, which is not in accord 
with the activation energy reported in this study. Over the 
temperature range of this present work the contribution to 
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the apparent diffusivity using the bulk value of diffusivity 
reported by Johnson [38] is negligible. 

If surface diffusion is assumed to be dominant, then a 
similar discrepancy exists between the reported activation 
energies for surface diffusion and those of the present work. 
If we consider the terrace ledge kink model of a solid-vapor 
interface developed by Kossel [59] , then the activated 
process for transport along a (111) FCC surface involves two 
s teps . 

The activation energy is considered to be the sum of two 
contributions, a change in enthalpy necessary for formation 
of the diffusing species and secondly, the change in enthalpy 
necessary for motion. The adatom in general carries the 
surface flux. It is the most weakly bonded atom to the 
surface, and on a FCC surface is characterized by 3 nearest 
neighbor bonds. An adatom is formed by an atom escaping from 
a kink site in a monatomic ledge (6 nearest neighbor bonds) 
to a site at the ledge (5 nearest neighbor bonds) and then 
to an adatom site. This involves the breaking of 3 bonds. 

The motion of an adatom involves the breaking of one bond. 

The activation energy therefore, can be looked upon as a 
process involving the breaking of 4 nearest neighbor bonds. 
The heat of sublimation of a FCC metal is the enthalpy change 
when a mole of atoms each atom breaking 6 bonds, escapes to 
the vapor state. Therefore, we can estimate the activation 
energy per mole of atoms as 2/3 AH su bii ma tion" For 2 old and 
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silver, the heat of sublimation is approximately 80 and 60 
kcal/mole respectively. This high value of activation 
energy is in accord with what has been reported for gold and 
silver where T/T ... is greater than 0.77. 

It has been suggested by Gjostein [51] that for studies 
at low temperatures (T/Tm < 0.77) the activation energy is 
approximately one-fourth the heat of sublimation. One 
study [60] of grain boundary grooving in gold reported an 
activation energy of 17 kcal/mole for T/Tm = 0.57. 

In order to be in accord with the presently accepted 
terrace ledge kink model, it must be postulated that the 
surface is covered with a large fraction of growth sites 
available for surface transport. The activation energy can 
then be assumed to involve the breaking of one bond. This 
model is for the most part unsubstantiated and sufficient 
experimental work has not been done in the present study to 
support or challenge this model. 

If we consider grain boundary diffusion to be the domi- 
nant mechanism, then we again are faced with interpreting the 
low value of activation energy experimentally observed (10,900 
cal/mole) . Hoffman and Turnbull [36] reported an activation 
of 20,200 cal/mole for a study of grain boundary diffusion of 
radioactive silver in a silver bi-crystal. If we just take 
into account the grain size of the film, the ratio of 
Dapp/DcB equals 0.01, 0.15 and 0.40 for temperatures of 200, 
300 and 416°C respectively. Let us assume the results of 
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Hoffman and Turnbull [36] as representative of a lower limit 
on the degree of disorder in a grain boundary. An upper 
limit would be to consider the grain boundary as in the 
liquid state. Activation energies for liquid self diffusion 
in silver have been reported [58] between 7 and 8 kcal/mole. 
The ratios of the apparent diffusivity to diffusion down 
"liquid" (highly disordered) grain boundaries are greater 
than 1 for all annealing temperatures. If we assumed that 
the grain boundaries between the grains in the thin film 
couples are highly disordered as the photographs of the film 
grain structure in Figure 25 and Figures 26 through 28 of 
the film surface suggest, then we would expect that the 
activated energy for transport along the grain boundaries 
to be between the upper and lower limit discussed above. 
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CONCLUSIONS AND RECOMMENDATIONS 
Conclus ions 

The study of lateral diffusion in thin film couples of 
gold and silver can be investigated using one of the non- 
steady state solutions to Fick's second law. The thin strip 
approximation, which in the limit of strip thickness = 0, 
is more commonly applied to a pair of semi- infinite solids 
can also be applied to a thin film of constant thickness 
where the initial profile is experimentally determined. 

The use of Zeibold's empirical conversion parameter, 
more commonly applied to the determination of composition 
in the bulk, is valid for thin films of the gold-silver 

o o 

system up to 2000 A in thickness. Up to 2000 A the thin 
film approximation that intensity is proportional to mass 
thickness is valid in the gold-silver system. 

Studies of lateral diffusion in Ag-Au couples exhibit 
surface morphology changes similar to that reported in bulk 
couples. The mechanism, operative in the thin film diffusion 
couple giving rise to the formation of a zone of porosity and 
protrusions on the silver and gold sides of the interface of 
the couple respectively, can be attributed to the conclusions 
reported by bulk investigators. The maintenance of continuity 
of the silver and breakup of the gold sides of the thin film 
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couple follows the predictions of Mullin's theory of grain 
boundary grooving. 

The use of incremental couples with small initial con- 
centration differences eliminates the formation of porosity 
and hillocks as experimentally observed in the bulk studies 
of Johnson. The maintenance of continuity in film thickness 
across the diffusion zone allows the determination of diffu- 
sivity as a function of temperature. Using the standard 
rate equation for an activated process the apparent diffusiv- 
ity was written as D = 1.9 x 10" 7 exp ■ cm 2 sec' 1 where 

W 4 

D 0 = 1.9 x 10 -7 and the activation energy is 10.9 kcal/mole. 
In addition a value of diffusivity was determined on the 
Ag-Au couple which had film separation by considering the 
porous zone as creating a new set of boundary conditions for 
the experimentally determined initial profile. The value of 
diffusivity was in accord with that observed in the incre- 
mental couples. 

A comparison of the possible transport mechanisms 
operative in the thin film couples between 200 and 416°C 
suggests that grain boundary diffusion is the dominant 
mechanism. In order to resolve the low value of activation 
energy reported in both the incremental couple and the Ag-Au 
couples an assumption that the grain boundaries were highly 
disordered must be made. 
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Recommendations 

In order to make a stronger case for the mechanism of 
grain boundary diffusion, incremental couples should be made 
with varying rates of deposition in order to vary the 
resultant grain size. Larger grains with a rising value of 
activation energy in accordance with the previously reported 
bulk values of grain boundary diffusion would lend more 
credence to the disorder assumption made in this study. 

A study of incremental couples with larger concentration 
differences could prove interesting. The magnitude of the 
initial concentration difference necessary to cause the 
formation of porosity could produce some interesting con- 
clusions as to the operative driving force for the formation 
of the porosity and hillocks. 

Because of the small amounts of material present in 
thin films and short annealing times, a study of lateral 
diffusion using the scanning electron microprobe could prove 
valuable. The small beam size would allow greater resolu- 
tion, in the lateral direction, and in the relationships of 
surface morphology to chemical composition, than in the 
electron microprobe. 

A followup study on the application of the empirical 
conversion parameter of Zeibold to other thin film systems 
could prove to be a valuable tool in film composition deter- 



minations . 
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The study of non-miscible systems using the experimental 
deposition technique of the present work should prove 
superior to the overlap technique used in previous work. 

The elimination of the overlap region removes the necessity 
for diffusion normal to the substrate to homogenize the 
overlap region and more closely approximates the initial 
conditions of the bulk sandwich couple. The method of 
deposition used in the formation of the incremental couples 
could be applied to the study of diffusion in non-miscible 
systems in single phase regions. 
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APPENDIX A 

A STUDY OF LATERAL DIFFUSION IN 
THIN FILM COUPLES OF BISMUTH AND TELLURIUM 

This author's initial investigation of lateral thin 
film diffusion took place in the bismuth- tellurium system. The 
separation between the razor blade and the substrate was 
achieved using a differential screw driving mechanism [9). 

The adjustment was done manually and resulted in regions of 
mixing on the order of 10" 5 meters. 

O 

All Bi-Te thin film couples were 1,000 A thick and 
deposited at room temperature. The anneals were carried out 
in two separate environments. Vycor encapsulated tubes at a' 
residual gas pressure of approximately 50 microns were used 
for anneals at 1, 3, 10, 30 and 96 hours. In addition, 
anneals in high vacuum at pressures less than 10' 7 torr were 
carried out for times of 1 , 3, 10, and 30 hours. 

The initial profile was determined on the electron 
microprobe, on an unannealed thin film couple. The region 
of mixing was 84 microns. From optical observation, and 
the electron microprobe intensity data, 71.5 and 12.5 
microns were attributed to the bismuth and tellurium overlap 
region respectively. The interface (position of razor blade) 
separated the two overlap regions. 
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Jn the series of anneals in the vycor tubes, the follow- 
ing observations were made: A peak in tellurium intensity 

was observed to move through the region of mixing towards 
the pure bismuth side of the couple into the tellurium over- 
lap region. The anneal at 10 hours signaled the completion 
of the shift of the tellurium "pulse," which now was entirely 
in the original tellurium overlap region. In the bismuth 
overlap region (tellurium side of original interface) a 
reduction of tellurium on the pure tellurium side of the 
couple was observed, followed with an increase in tellurium 
over the entire bismuth side of the couple. The pure 

O 

tellurium film thickness originally at 1,000 A was reduced 
to 200 A. 

The observations of the couples annealed in high vacuum 
were the same for 1, 3 and 10 hours. However, the 30 hours 
anneal exhibited no significant decrease in tellurium 
thickness on the pure tellurium side of the couple or 
subsequent increase in tellurium intensity over the bismuth 
side of the couple. 

The bismuth intensity profile remained immobile for all 
annealing times in both environments. X-ray fluorescent 
analysis using Cu Ka radiation confirmed the presence of the 
intermetallic compound Bi 2 Te 3 in the region of mixing for 
annealing times greater than or equal to 3 hours. After 10 
hours no further growth of the intermetallic compound was 



observed . 
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The mean free path of ambient gas atoms at 50 microns 
is 0.5 cm. The rate of free evaporation of tellurium and 
bismuth at 200°C is 0.2 and 2 x 10' 5 monolayers/sec., 
respectively. Therefore, mass transport by the vapor phase 
is a probable mechanism in the couples encapsulated in the 
vycor tubing. This would explain the increase and decrease 
in counts of tellurium over the bismuth and tellurium sides 
of the couple for the 30 and 96 hour anneals in the vycor 
capsules . 

The formation of the intermetallic compound occurs in 
both environments, ruling out a vapor transport mechanism, 
especially unlikely in high vacuum at 10' 7 torr where the 
mean free path is 450 meters. 

The relative immobility of the bismuth profile, and the 
moving tellurium "pulse" for the anneals at 1, 3, and 10 
hours in both environments led the author to conclude that 
the tellurium was the only active (within the time frame of 
the author's study) species in the mass transport process. 

Let us consider the region of mixing to be made up of a 
series of thin strips side by side in the lateral direction. 
Because of the simultaneous deposition process, each strip 
will be homogeneous in composition in the direction normal to 
the substrate. The variation in strip composition in the 
lateral direction should be from pure tellurium to pure 
bismuth across the region of mixing starting from the pure 
tellurium side of the region. Therefore, one strip in 
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original region of mixing should be of the stoichiometric 
composition of the intermetallic compound Bi 2 Te 3 . If we 
arbitrarily assign the mass ratio of bismuth to tellurium 
(m Bi /m Te ) in this strip of stoichiometric Bi 2 Te 3 equal to 1, 
then strips on either side would have a mass ratio < 1 and 
> 1 towards the pure tellurium and bismuth side respectively. 
Mass transport of tellurium in the direction of the pure 
bismuth side of the couple would increase and decrease the 
respective mass ratios towards 1, on the tellurium and 
bismuth sides of the initially stoichiometric strip. The 
region of mixing would now be a series of strips of Bi 2 Te 3 . 

However, after the 10 hour anneal which signals the 
conclusions of the mass transport process, the region of 
mixing, and now zone of intermetallic compound, appears to 
remain stationary for further times. If the compound forma- 
tion was diffusion controlled, a parabolic growth of this 
zone with respect to time, would be expected. Since some 
other time law is governing its process, it would appear to 
be an interface controlled process. 

This, however, creates an anomaly in that mass transport 
was observed across this small strip of originally formed 
intermetallic which certainly had some sort of interface 
associated with it and yet the interfaces of the strips on 
each side of the region of mixing limit the intermetallic ' s 
growth to some small time power law. 
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A possible explanation for this behavior is that surface 
diffusion is the dominant mechanism. Instead of tellurium 
atoms moving across the interface, consider mass transport 
across the surface of this small strip of initially stoi- 
chiometric Bi 2 Te3. This causes tellurium accumulation on 
the bismuth side and mass depletion on the tellurium rich 
side of the strip. Followed by diffusion into the film 
perpendicular to the substrate, the stoichiometric strip 
will grow laterally. However, as the strip grows in lateral 

i 

extent, the distance of tellurium surface transport would 
increase. If the rate of transport across the surface 
decreases as the extent of the strip grows, as is experi- 
mentally observed when the extent of tellurium transport is 
limited to the region of mixing, then the rate of formation 
of the intermetallic is fast initially falling at an unknown 
decreasing rate. 

This strip as it grows would also provide a barrier to 
the tellurium moving away from the tellurium rich side of 
the region mixing, which at the same time, through tellurium 
depletion, is creating strips of intermetallic, even though 
very thin in thickness normal to the substrate. 

With the intermetallic extending throughout the region 
of mixing, and now due to its extent in lateral dimension, 
and its respective interfaces, we would not expect and do 
not observe any growth with respect to time. 
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In summary, we can conclude that thin film couples of 
bismuth and tellurium form the intermetallic compound in the 
region of mixing and further growth of the intermetallic 
compound was not observed. Therefore, the intermetallic 
compound acts as a diffusion barrier and transport of bismuth 
and tellurium through the intermetallic is interface con- 
trolled following some small time power law which was not 
determined in this work. Due to the vapor pressure of 
tellurium, it is recommended that any further studies in 
the bismuth- tellurium system be done in a high vacuum condi- 
tion and/or temperatures below 200°C. 
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APPENDIX B 

EXPERIMENTAL RELATIONSHIPS USED FOR DEPOSITION OF 
ALLOY FILMS AND THIN FILM COUPLES 

Figure 38 shows the orientation relationship between 
the source, substrate and crystal monitor. The equations 
for evaporation from a small area (source) onto a parallel 
plane receiver had been previously reported [9] . The change 
in frequency (frequency shift) of the 5 MC crystal is related 
to the mass deposited on the crystal by equation (25) of the 
present work. The 6 variable represents the correction param- 
eter that adjusts for the orientation relationships of 
Figure 38. Therefore, the t in equation (25) represents the 
thickness of the film. 8 consists of the following variables 



R! = the substrate to source distance 
R 2 = the crystal to source distance 
0 = angle between the source normal and the crystal 
monitor . 

The thickness variation across the film is 



^ 2 (R 2 J cos 2 0 




(30) 



where j is a constant of the crystal 



to 



1 



( 31 ) 
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Figure 38 

Orientation Relationships between the Source, 



Substrate and Crystal Monitor 
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where t = thickness of film directly above source 

t 0 = thickness of film at distance r from point on 
substrate directly above source. 

Using equation (31) the maximum thickness variation for a 
1,500 X film using the values of r = 2 cm, R = 18 cm, used 

O 

in the present study, is 38 A, which is within the error of 
thickness measurement of the Dektac. 

Alloy deposition was carried out by 1) assuming that 
the film density was that of the reported bulk value and 
2) that the total thickness of the alloy film, t, was equal 
to the sum of the thickness of the two pure films. There- 
fore, let 



t 



'Au 




(32) 



Let desired composition of film be X weight! Ag , (1-X) 
weight! Au. 



X weight! Ag = 



mas f A 8 X 100 

mass Ag + mass Au 

p Ag t Ag A 

p Ag *"Ag A + Pj\u *"Au A 



(33) 



where A = substrate area on which gold and silver is 
deposited . 



Substituting 
for t . 

Ag 



equation (32) into equation (33) 
tp A u( w S Ag) 

t Ag ' P Au (Vn + p Ag (W% AuJ 



and solving 

(34) 



t 



Au 




and 



t 



(35) 
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APPENDIX C 



THIS IS R PROGRAM T0 DETERMINE THE EMPIRICAL ZEIBOLD 
CORRECTION PARAMETER FROM THE THIN FILM STANDARDS OF 
THE SILVER- COLD SYSTEM 

»1 LINEAR LEAST SQUARES FIT IS DONE FOR EACH ALLOY 
COMPOSITION MADE AS A FUNCTION OF FILM THICKNESS. 

FROM THE RESULTANT EQUATIONS DATA IS EXTRAPOLATED FOR 
FILM THICKNESSES OF 500, 1000.. 1500, AND 2U00 ANGSTOMS. 

THE VALUES OF K AND C ARE DETERMINED AND A cl NEAR 
LEAST SQUARES FIT IS DONE FOR A PLOT OF 1-K/K VRS. 

I -C/C, FOR EACH FILM THICKNESS. 

FROM THE SLOPE THE ZEIBOLD CORRECTION PARAMETER IS 
DETERMINED. 

COMPOSITION FOR ALL THE DATA IS REPORTED IN ATOMIC 
PERCENT SILVER 

K IS THE BACKGROUND CORRECTED INTENSITY/BACKGROUND 
CORRECTED PURE INTENSITY 

C IS HEIGHT FRACTION COMPOSITION 

BKGD1 IS THI BACKGROUND OF THE PURE ELEMENT 

BKGD2 IS THE BACKGROUND OF THE ALLOY 

I NT IS THE TOTAL INTENSITY FROM THE MICROPROBE 

THK IS THE FILM THICKNESS 



DIMENSION INTCIO), THK 10 > , AOCICO, HI <10>, I PURE <5, 5), I ALLOY <10, 10>, K 

j a- io, io>, cci co, xcio>, ycig, io> 

REAL KFl, 1 PURE, I ALLOY, I NT 
B"0. 0 
. 1=1 

1 IFCJ EG 1 ' HRITECG, 112;' 

IF<J EC! 2 WPITFC6, 1.13.' 

D01.0 L=l, 7 

N=4 

IF CL. GT. :1> N=5 

READ <5, 100> COMP, (I NT' I > 1=1, N> 

READ <5, 101 > f THK CO, 1=1, N> 

IF'.N EQ. 4 AND. J. EQ. 1) GO TO 20 
IF. N EQ. 4 AND J. EQ 2> GO TO 22 
GO TO 21 
ij HRITE' G, 1 14 > 

GO TO 23 

2 HRITECG. H5> 

GO TO 22 

1 HRITECG, 102> COMP 

3 HRITECG, :t TO - 
UR I TECS, 103 > I NT < I >, THK' I >, 1=1, NT 
SUM1=0. 0 

3UM2=0. 0 
3UM3=G 0 
,UM4=0 0 
3UM5=0. 0 
DO 1 1 1=1, N 
:Y=IrJ R 0*TEIM I> 
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V V-- 1 NT < I ) + INT a> 

XX«THK<n*THK<I> 

SUM1=THK<I>+SUM1 
SUMS* I NT ( I > +SUM2 
SUM3=XV+SUM3 
SUM4=XX+5UM4 

11 SUM5=YV+SUM5 

RN- < SUM2+SUM4 > - < SUM1+SUM3 > 

RD= ( N+SUM4 > - •! SUN1+SUM1 > 

R1N= < N+SUM3 > - < SUM1+SUM2 > 

R0<L)= s RH/RD 
R1<L>=R1H/RD 
MRITE<6, 111) 

S YX1=SUM5- < R0 C L ) +SUM2 > 

SVX2= SYX1- < R1 C L. ) *5UM3 ) ) /to 
SYX=SQRT<SYX2) 

S I G1=SUM5/N 
S I G2= < SUM2/N > + < SUM2/N > 

S V2=S I Gl-S I G2 
SV-SQRT < SV2 ) 

CORR=SQRT (1. G-CSYX2/SY2) ) 

NR 17 E ( 6, 107 R0 < L > , fil < L > , SYX, SV, C0RR 
IF a EQ. i> GG TO 10 
J FCL GT. 2) GO TO 12 
THICK -500. 

D013 M=l.< 4 

BKGDl=R0 1 ) + (Hl(l) +THICK) 

I PURE <2. M)=<fil<L)*THICK> + R9<L) - BKGD1 

13 TH I CK=TH I CK+50Q. 

GO TO 10 

12 THICK=500. 

IF <L.. EQ. 3> GO TO 10 
D014 M=l. 4 

BK=RG<3> + <R1(3)*THICK> 

IF<J. EQ. 2) GO TO 40 

BKGD2-BK+ < 1. -COMP) + BKGDl+COMF 

GO TO 41 

40 9 K G D 2 = BK +C O M P + BKGD1*<1. -COMP) 

41 I RU. O V ( L.- M ) = ( fll < L ) +T H I CK > + RO(L) - BKGD2 
KRCL. M> = IRLL0V<L, M>/IPURE<2, M) 

T H I CK=TH I CK+500. 

14 V < L, M > = ( 1 . 0-Kfl < L, M > ) /KR < L, M > 

Cl “1 078 8* COMP 

C2-1. 97+ a. “COMP ) 

C < L ‘>“01/01:1 f C2) 

I F J EQ. 2) C(L.)=1. O-CCL) 

X<L> = <1. I>C(L.))/C<L) 

10 CONTINUE 
THICK =500. 

RLL'1-0. 0 

RLL2=Q. 0 

hLL.3=0. 0 

RLL4=0. 0 

DO 15 M=l, 4 

WRITE (6, 105;:- THICK 

TH I CK=TH I CK+500. 

WRITE (b, 1.04; 

D016 L=4, 7 

16 NR I T E ( 6, 106 > V ( L. M > - X ( L > 

SUM 1=0. 0 
SUM2=0. 0 
SUM 3=0. 0 
SUM4=0. 0 
D018 L=4.< 7 
XX»X<L)*X<L) 
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VV=V'L, M)*Y<L, M) 

XY=X(L)*Y<L, M) 

SUMDXX+SUM1 
SUM2=VV+SUM2 
SUM2=XY+SUM3 
18 SUM4=Y<L, MHSUM4 
ALL:l=RLLl+SUMi 
ALL2=ALL2+SUM2 
ALL3 =hLL3+SUM3 
ALL4=ALL4+SUM4 
SL.0PE=SUM3XSUM1 
SVX2= < SUMS- < SL0PE*SUM2 ) > /4. 

SYX=SGRT (SYX2) 

SIGl=SUH2/4. 

SIG2=(SUM4/4. )*<SUM4/4. ) 

SV2=SIG1-SIG2 
SV=SQF:T (SY2) 

C0RR=SQRT < 1. j- C SVX2/SV2 ) > 

15 WRITERS, 10?) B, SLOPE, SYX, SV.. C0RR 
5L0PE=ALL2/ALL1 
S VX2= < ALL2- < SL0F'E*ALL2 »/lS. 

SVX=SQRT < SVX2 > 

SIGDALL2V16. 

SIG2=<ALL4/1S. >*<ALi_4/16. > 

SV2=SIG1~SIG2 

5Y=SQRT(£Y2> 

C0RR=SQRT < 1. G- < SVX2/SV2 ) ) 

WRITE <6, 108 > 

WRITE <6/ 107 > B, SLOPE, SYX, SV, C0RR 

IFCJ EQ. 2) STOP 

J=2 

GO TO 1 

100 FORMAT <F5. 1, 5X, 5F10. 5) 

101 FORMAT <5F 10. 5) 

102 FORMAT <1H0> ,’///> 20X, ' (ATOMIC FRACTION SILVER)' , X, 24X, ' COMPOSITION 
1= ', F5. 2) 

103 FORMAT < 2 < SOX, F13. 6» 

104 FORMAT (1H0, 17X' 1-K/K' . SOX, ' 1-C/C' ) 

105 FORMAT d.HO, ////> . SOX, ' FILM THICKNESS = ' , FIG 2) 

IOC FORMAT <2d5X, FID. 6)) 

107 FORMAT < 1H0, 20X, ' INTERCEPT (AO) = ' , 26X, F13. 5, /V, 21X, ' SLOPE <AD = 

.1/, 30X, FI 3. 5, //• 21X, 'STANDARD ERROR OF THE ESTIMATE = ' , 10X, F13. 5, / 
2/j 21X, 'STANDARD DEVIATION OF V = ', 17X, F13. 5, //, 21X, ' COEFFICIENT 0 
3F CUk'REl H I ION = , 14 X, FI 5) 

108 FORMAT OHO, ////> 15X, 'LINEAR LEAST SQUARES FIT FOR ALL DATA') 

109 FORMAT (II ID 

110 FORMAT dHO, 22X, ' INTENSITY', 22X, 'FILM THICKNESS' ) 

111 FORMAT ClHO, 15X, 'EQUATION OF LINE IS INTENSITY = AO + R1 X T 

THICKNESS ') 

112 FORMAT (Ihl, 15X, ' CALIBRATION DATA TAKEN OFF THE SILVER SPECTROMETER 

1' ) 

113 FORMAT ( 1H1 , 15X, 'CALIBRATION DATA TAKEN OFF THE GOLD SPECTROMETER' ) 

114 FORMRTdHO, 15X, 'BACKGROUND INTENSITY OF PURE SILVER FILMS' ) 

1:15 FORMRTdHO, 15X, 'BACKGROUND INTENSITY OF PURE GOLD FILMS') 

END 



APPENDIX D 

COMPUTER PLOTTING PROGRAM WHICH PLOTS CONCENTRATION 
VERSUS RELATIVE DISTANCE FROM 



MICROPROBE INTENSITY DATA OF THIN FILM COUPLES 
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APPENDIX D 



THIS IS H PI Ain I NG PROGRAM WHICH FIRST DETERMINES THE 
CONCLN (RATION r ROM THE IN TENS I TV OUT OF THE NfCROPROBE 
HMD PLOTS CON CENTRA! ION VERSUS RF LA f I VE DISTANCE 

BKGD1 IS THE BACKGROUND CORRECTION FOR THE PURE FILM 

BKGL'2 IS THE BACKGROUND CORRECT JON FOR THE ALLOY FILM 

no HR is I he sfjoold parameier determined from INDEPENDENT 
MEASUREMf NT S OF THIN FILM S f RNDERDS 



PURE I IS IMF Rl ll-l INTENSITY OF THE PURE ELEMENT MPhSURED 
FROM THE; MIiROPROBE SPEC I ROME TER 

RL-.V < I > is THE RUN INTENSITY OF THE ALLOY THIN FILM MEASURED 
DIRECTLY FROM THE MI CROPROBE SPF C TROMET E' R 



HL.DE N IS THE ALLOY DENS I TV FOR THE SILVER - GOLD SYSTEM 
COMPUTED ON I HE BASIS OF HN ASSUMED LINEAR RELATIONSHIP 
THE ALLOY DENSITY AND THE ATOMIC FRACTION OF THE PHASE 



AND I 
BETNEhN 
SYSTEM 



D I MENS ION Hi. V 5 AQ > , CONC 500 > ■ D I ST C 500 > , CONCi ( 5UU > .. D I ST 1 < 500 > , C6NC2 
1- 500 >, D I S I 2 v jUO j .. I MAGE <150 A , CTESKIOO), PTESTOiOG), NSCALE<5> 

REAL. K, Km. KAlT 
INTEGER PUN, Z. V, G 

nscalk<. i > - j. 

r JSCAI T 1 2 “0 
N SUM I F 2 

NSCAL F < 4 > 0 
r JSCAI F* 5 > -• 1 
PI =2. 1.4 ISO 
PX -C L -05 

• Ml I PLOT 1 <NSCAI. E f 5, 10, 1.1 . :1 0> 

1 REAL 1 < 5, I i in .■ RUr J • N, A BAR, BKGD1 , BKGD2, PURE I , T I ME ■ I L: MR 
IT 1. 

I.ji" 0 

DO in, n 

READ <5, .1.02 > ML V < O 

2 PIS T < I ) “U B-KFLOlira - l > 

L - PURE J Bl i.iL'1. 

1-1 

9 D“ HLV< I ; 'BKGD2 

IF l F. 0. CO GO TO 3 
K- D/E 

\ H = <1 . U -K '/(LaABAR) 

I MR- L 0 A: A 
C=l. G/KAfl 
GO TO 4 

2 f ; IJ IJ 

C-0. 0 
GO TO 4 

4 HG" C/i. 0Zo8 
flU- <1 -C :■ /l. 9P 
A T MF R ; HU, ’ ■ RG i FlU 1 
ML PLNXL9. 415 - <3. 843 + AT MFR • 

L Of JC < I >=ALPF.N+C 

I F >, i.;ONC < I > . G T 1 56. OR. CONC<I •. I T 0. 56 > GO TO 20 
Ci ”0 + L 

C'TEST • G>- CONC< I 1 

dtp* t < g ’ Disra :• 

20 IFa-ll:l:' 5, 6. 6 

f. I Hi I PL0T2< IMAGE, 1500, SO 0, 0. 0, 2. 5, 0. 0> 



in f 
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< m i. pi pi i • < iih-, rvisr, cgnc, i> 

14: 1. IT •: A, I in :> 

f’l Mi ,p,'H CPNCf.N'f Rfl f I ON IN GROMS RG /CO 

MF-' uf;. p, i iT.r.. 

G0 10 J5 

5 IF ( I — N '• 12, 6, 15 
12 I=Itl 
GO TO 9 
15 SUM 1-0. 0 
SUM2=0. 0 
SUM 2=0 0 
SUM4=0. 0 
SUM5=iJ 0 
D021 1=1, G 
XV=C'T EST < I >*OTEST< I 
XYV=C'I ES Tc I > f CTEST < I > 

XX=D1ESTC I r ST< I •' 

SUM1=D TE ST ( I :> \ :;i.JMl 
SIJM2=C FESra> n.;i IM2 
SUM2=XV+GUM2 
SUM4==XX+SUM4 
21 SUM5=XYY4SUM5 

AN= C SUM2 I SI IM4 > - SUMl+SUM.l: ) 

HD= < FLOAT G > A SIJM4 > -• SUM1 + SUM 1 > 

Hi N~ < FL0R T r. G /■ •+: SUM 3: ) - < SUM1 +SUI12 > 

RO=flN/AD 

Al=fllN/AD 

SVX1-SUM5- < HO f SUM2 > 

syx 2 =csv::i-< hh sumo vflgat cg> 

5VX=SQRT<SYX2> 

SIG1=SUM5/F1.0RT(G> 

SIG2=”<SUM2/F I ORT < G > > + <SUM2/FL0AT XG> > 

S V 2 = S I G 1 - 5 1 G 2 
SV=SQRT<SY2> 

CORR~SOR f ( 1. 0“ C SYX2XSY2 > > 

NR I TEC 6, 1 1 1 > RUN, T I ME, TEMP 
NR I TEC 6, 1 12> Hi, HO. SYX, SV, CORR 
GO TO 1 

100 FORM H K IT, 13, 4X, SF10. 5> 

101 FORMATCiHl , IX, ' RUN NO ", 5X, RBRR = ", 8X, 'BKGD1 = " , 7X, " BKGD2 = ',? 
IX, -'PURE I - ",?X, "N - " , 5X, •' T I ME = " , 7X, " TEMP = ' > 

102 FORMRT C IMG, 15, 4X, F8. 5, 2C5X, F10. 1 >, 7X, 12, 5X, FC. 1, 8X, F6. i> 

1 1 1 FORMR T Ti l 5 ..- 

1 04 FOR MR “I < 1 HO, 4X, RLLO'T I NTfc NS 1 T V" , 5X, " KAG , 9X, " NT FRRCT. " , ?X, " ATOM 
1FRRCT. " . 5X, " RLLOV DENSITY*- , 2X, " CONCENTRATION" > 

105 FORMRT '• 6 CSX, F10 5>) 

inf f-ORMRTMHl, 5X. "DIFFUSION CRFFF Tr TENT FOUR IS " . 1 PE12 4> 

•i i tv r Mpi !RT' /t: .' r i o v. 

I U1 • I I'll ! I 1111' '.I,!, UT'T J. ! 11 ■ I I It'll • t , ti-'-i, (.I! 1 I .' 

in; i hfiu ; I •: in . * •' c : • rmuiivi i.'ImiiNi > M 1 1 y: nr . • * •* 



11 U I ULl 'll ! T \ J. 1 II. 

111 FORMRT < 1H 
1 ANNEAL = " 

112 FORMRT <1H 
II ON AXIS = 

20X, FIT. 

3TCIENT Of COF’PEI HTluN = 

1 12 FORMRT V 2 OX, ? < 10X, F8. 2 > > 
114 FORMRT Cl HI, 20X, "G = ", l. 



, Ft. 1, 10X, " TEMPERATURE OF 



20 X, " RUN = ", 13. 10X, - TIME 

re. i> 

SOX, SLOPE = " , 35X, FIT. 5, /, 21X, " INTERCEPT OF CONCENTRRT 
. J1X.F13. 5, /, 21X, " STANDARD ERROR OF THE ESTIMATE = M 
, 21 X, STANDARD DEVI AT ION OF V = " . 17X, F13. 5, /, 21X, " C6EFF 
14 X, FIT. 5) 



END 



APPENDIX E 

COMPUTER PROGRAM WHICH NUMERICALLY SOLVES EXPERIMENTAL DATA 
OF THE INITIAL CONCENTRATION PROFILE FOR GIVEN VALUE 



OF D AND TIME 
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APPENDIX E 

c 

c 

c mis is tmc numerical test program which operates on the 

C EXPERIMENT AL DR T A 0F RUN 22^ USING THE THIN FILM S0LUTI0N 

C 

C RFJRR IS THE 2EIB0LD PflRRME T FiP DETERMINED FR0M INDFPENDRNT 
C MERSUREMEN TS OF THIN FILM STRNDERDS 

C 

C PURE I IS THE PAW INTENSITY OF THE PURE ELEMENT MEASURED 

C FROM THE MICPOPROBE SPECTROMETER 

RLV ( I > IS THE RRW INTENSITY OF THE RLL0V THIN FILM MEASURED 
C DIRECTLY FROM THE MJCROPROPE SPECTROMETER 

C 

C BKGD1 IS THE BACKGROUND CORRECTION FOR THE PURE FILM 



BKGD2 IS THE BACKGROUND CORRECTION FOR THE ALLOY FILM 

ALDEN IS THE ALLOY DENSITY FOR THE SILVER - GOLD SYSTEM AND IS 
COMPUTED ON THE BASIS OF AN ASSUMED LINEAR RELATIONSHIP BETWEEN 
I HE ALLOY DENSITY AND THE ATOMIC FRACTION OF THE PHASE SYSTEM 



ATMFR IS THE ATOMIC FRACTION 



DIMENSION ALY0500), DIST<500>, C9NC<500>, TOCONC*: 500;a CODELTCIOCO, TCG 
INC < 100), IMAGE* 1500), DISTK 1.00), TEST *100), NSCRLE»'5> 

REAL K, KA, FAR 
INTEGER W, RUN. 7, Y, 0 
Y-0 

PI -3. 14159 
DX=S. E- 05 
NSC ALE 1 > = 1 
NSCHl.E«:'2>--0 
NSCALE 
N SCALE < 4 )=0 
f -I SCALE 05 ,rr 2 

1 RL HD CEL TOO RIJN, N, REAR, BKGD1. . 8KGD2, PURE I, TIME, TEMP 
API TE-S, 10 1 > 

TIP J 1 EOS, 102) RUN. RUHR, FJf GDI , BKGD2, PURE I, N, TIME, TEMP 
r f. ML 48. 0 
DO 2 I -i. N 
RL AD 1 - 5, 10?) Al Yd) 

2 D1SKD-0 o+FLOAT < I -l) 

E-PUF'EI -LI GD L 

t t 

TIP J TEC G, J u4 > 

9 D-ALY'.. I )- BKGD2 

IF CD . LF 0 0) GO TO 3 
E>D/E 

If <K. GE. 1 CO GO TO 14 
KA - Cl OdOAK+ABRP) 

KRA=i O+KR 
C=l. O/KFIA 
GO TO 4 

3 L-0. 0 

c=o. 0 

GO TO 4 
14 K-l 0 
i> 1. 0 

4 AG= C .01 . 07A8 
MU -a. -0/1 37 

fl I MFR- AG/ * AG i- AU ) 

ALDFN=t9 415 - 03. 843+AT MFR > 



CONC<I)=RLntN'iC 

WRITEC6, J 05 > RLV< I >.. K, C, flTMFR, fiLDEN- C0NC<I ) 

5 1 0CONC < I > ~C0NC < I 

6 IF 'A. EQ. N> G0 TO 8 
I-Ifl. 

GO TO 0 
8 M N i 
DO 10 I I. iT 

CiiDELTc I :• ~ < T CiCONC ( I > +T QCfTNP < I+l)}/2 
10 CONTINUE 
99 DIFC=5. F- 11 
98 W- 0 

RDT2~2. -KSQRT C D I F 03600 fT I ME ) 

C' T4-0IFC + 1 IME+O. 144E05 

ROOT L> l “SOFT < p I I F C+ 2600. +T I ri E > 

DO 12 Z-i, M 
SUM-0. 0 
D013 L=l. M 

DEI T - < FLO.TT ( Z > - FLOAT < L ) ) +DX 
DEL TSQ -DEL T +DEL..T 
0 1 V~~ < DEL TSQ/DT 4 ) 

E=EXP ( D I V > 

PROD-E+CODF.LT < L > +DX 
SUM-SUM-* PROD 
13 CONTINUE 

FUN- < FL OR r < 2- 1 ' +DX > /RDT2 
rC0Nr:l-<2 12/2 )*ERFC<FUN> 

T C 0NC2- SUM/ < 2. +ROOTD T > 

T C ONC < Z > = T CON Cl f TC0NC2 

IF<TC0NC<7*:« GT. 1 56 OR. TCQNC(Z). LT. 0 56) GO TO 12 
W-W+l 

T EST ' W > ~T CONC < Z > 

D I STIC U^Di ST 
12 CONTINUE 

HP I TEC 6/ 106) PIFC 
WRITE <6, LJ 0/ TIME 
I" HI L FT on ■' NSrm F, 5, 10, 10, 10 > 

cull plot 2 < image, 1500, so. 0, 0. 0, 2. 5, 0 o> 
cull plot2*. ih*, dist, tconc, z*> 

COLL PL0T4<42, 42H CONCENTRATION IN GROMS HG /CO 

WRITE< 6. 109) 

V" V+l 

IFCV EQ. 1) TIME -12. 0 
IF<V FO 2) ST OF* 

DirC-5. 9E-11 
CiO TO 98 

100 FORMAT <15* 13, IX, CF 10. 5) 

101 F ORMAT < 1H1, IX, " RUN NO. ' , 5X, " RBRP = ' , 8X, " BKGD1 = ' , 7X, ' BKGD2 = ',? 
IX, ' PURE! - " , 7X. " N = '* , 5X, "TIME = " , 7X, " TEMP = " > 

102 FORMAT <1H0, 2X, 15, 4X, F3 5, 3<5X, F10. i>, 8X, 13, 7X, F6. 1, 8X, Fb. 1> 

103 FORMAT*:; FI 3. 5> 

104 FORMAT ■' 1H0, 4X, ALLOY INTENSI TV" , 5X, "KAG", 9X- " WT. FRACT. ", 7X, "ATOM 
1 TRACT . " , 2X, " ALLOV DENS I TV" , 2X, " CONCENT RAT I ON > 

105 FORMAT <6<5X, F10 5>> 

106 FORMAT am 5X, "DIFFUSION COEFFICIENT EQUALS 1PE12. 4) 

107 F 0RMR1 * 5X, F10. 5> 

108 FORMATTING* 5X, ' CONCENT RAT ION' /5X, " GRAMS/CC " > 

109 FORMAT < 1H0, 55X, "RELATIVE DISTANCE (MICRONS)-* ) 

110 FORMAT ( 1H0, r 55X, • TIME = ',F4 ; • 

111 FORMAT ( 10X. F10 5, 20X, FLO 5* 

112 FORMAT vlHO, 5X, " TEST CONCENT Ri! f l ON' , 5X, " RELATIVE DISTANCE (MICRONS) 
1" > 

11" FORMAT a HO, 20X, ' SLOPE = ", 25X, F13. 5, //, 21X, ' INTERCEPT OF CONCENTRA 
IT I ON AXIS = ’, 09X, FI 3 . 5, //, 21 X, "STANDARD ERROR OF THE ESTIMATE = ' 
2, 10X, F13. 5. //, 2 IX, " STANDARD DEVIATION OF V - " , 17X, F13. 5, //, 2 IX, X C 



APPENDIX F 



STATISTICAL RELATIONSHIPS USED IN THE PRESENT WORK 
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APPENDIX F 

STATISTICAL RELATIONSHIPS USED IN THE PRESENT WORK 



Consider curve of the type Y = aiX. Assume that for 
the data point (Yi, Xi) each Yi is an observation from a 
normal population with a mean aj Xi and constant variance 
o 2 . Assuming that measurement of Y locates it within a 
small interval of extent AY, the probability of observing 
the value Yi is from the normal distribution: 

N 

S = l (Yi - a x Xi) 2 (36) 

i = l 

3S 

For S to be a minimum — = 0 

3a i 



as 

aa x 



2EYX + 2a ! EX 2 = 0 



Therefore , 



EXY 
ai " EX 2 



(37) 



The variance about the regression line is estimated by the 
square of the standard error of the estimate 



2 _ E (Y - a t X) 2 

i X N 



EY 2 - ajEYX 
N 



( 38 ) 



where N = the number of data points. 
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Lines drawn parallel to the regression line at vertical 
distance (Y) of + Sy^., + 2S YY and + 3S YY will enclose 68%, 



'YX 



YX 



95% and 99% of the data points. The value of Sy-^ gives a 
measure of the dispersion about the mean. 

The coefficient of correlation is defined as 



r = 



1 - 



s 2 ' 
i X 



S 2 
Y 



(39) 



where 



2 _ £Y 2 



= 

b Y 



N 



IY) 2 

~n 



The coefficient of correlation is a confidence test on 
the linearity of the data. For example, if the coefficient 
of correlation equals 0.658 and 0.708 we can say with 98 and 
99% certainty that there is linear correlation of the data. 

If we consider a curve of the type Y = a 0 + aiX, then 
following the procedures above we can derive the following 
relationships 



(IY) (IX 2 ) - (IXI(IXY) 
NIX 2 - (IX) 2 



(40) 



NIXY - (IX) (IY) 
&1 NIX2 - (IX) 2 



(41) 



S 2 , the square of the standard error of the estimate equals 

I A 



■2 _ IY 2 - a 0 IY - ailXY 



YX 



N 



(42) 
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The variance of slope, V 



i 




( 43 ) 



The confidence band of the slope, say equal to b , is a fan 
shaped area converging on the mean value of X and Y with 
slopes b + Vt, where t is the probability level determined 
from an F test on the coefficient of correlation (equation 
39 ) . 
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